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Advances in the fabrication and characterization of nanoscale systems now allow for a deeper 
understanding of one of the most basic issues in science and technology: the flow of heat at the 
microscopic level. In this Colloquium we survey recent advances and present understanding of 
physical mechanisms of energy transport in nanostructures, focusing mainly on molecular junctions 
and atomic wires. We examine basic issues such as thermal conductivity, thermoelectricity, local 
temperature and heating, and the relation between heat current density and temperature gradient 
- known as Fourier's law. We critically report on both theoretical and experimental progress in 
each of these issues, and discuss future research opportunities in the field. 
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I. INTRODUCTION 

Understanding how heat is carried, distributed, stored 
and converted in various systems has occupied the minds 
of many scholars for centuries. Recently, the problem has 
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garnered even more attention and has grown considerably 
in importance. This is not due only to purely academic 
reasons: its practical impact in society has been recog- 
nized as one of the most critical programs for the devel- 
opment of the neces sary resources to sustain the future 
welfare of mankind 

In conjunction with these motivations, research seems 
to suggest that nanoscale systems (such as carbon- 
based nanostructures, organic molecules, etc.) may 
be good candidates for such technological advances. 
For instance, the flow of heat in nanosc ale systems 
may be ha r nessed yia therm o electr ic effects (Bel],[2Q0F 



iMaiumdarl . 12004 iRodgersl 120081) to generate heat- 
voltage converters, which (if their efficiency can be im- 
proved) may have real impact on global energy consump- 
tion. Other interes ting applications , such as nanoscale 
local refrigerators (|Shakouril. 120061). thermal t r ansis- 
tors (iFranceschi and Mingol. 120071: iGiazotto et all 120061: 



iLi et all 120061: lLo et »/]. 120081: ISaira et all 120071) ther- 



mal rectifiers iLi et a/1 l2004allbl: ISeeal and Nitzanl. [2005 



Terraneo et all l2002t IWu and Li 120071: lYang et al 
2009) and nanoscale radiation detectors (jGiazotto et al. 
2006) and e v en therma l memory and logic gates 
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(jWang andU I2007L l2008( ) add to the importance and 
interest of this research field. 

In spite of the recent advances, this research program 
still presents quite a few challenges related to the intrinsic 
non-equilibrium nature of the problem. In the presence 
of a heat current, quite generally, both electrons and ions 
may be very far from their equilibrium state. In addition, 
they are in interaction with each other and, at the same 
time, in dynamical interaction with one or more environ- 
ments. 

To complicate matters, heat flow is in many ways (as 
we will discuss in detail in the following sections) funda- 
mentally different from charge flow. Therefore, many of 
the theoretical tools which are used to describe charge 
transport cannot be straightforwardly and uncritically 
extended to the study of heat transport. From an experi- 
mental perspective, studying energy flow at the nanoscale 
is in several ways more challenging than studying charge 
transport, one reason being that no simple device analo- 
gous to an "ammeter" is at hand to measure energy cur- 
rents. Furthermore, the scale of achievable thermal con- 
ductivities is generally much smaller than that of electri- 
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cal conductivities ( Maiumd"arl 120041) . Consequently, one 
has to necessarily introduce models by which the thermal 
conductance can be deduced from measurable quantities 
such as charge current, voltage and temperature. In ad- 
dition, measurement schemes with macroscopic probes 
are necessarily used so that the channeling of heat only 
across the junction is difficult to achieve. 



(a) Quantum point contact 







J~10-10"nm 




~10 2 -10 3 nm 





(b) Suspended nanotube 




(c) Molecular junction 




FIG. 1 (Color online) Schematic representation of the dif- 
ferent systems we consider in this Colloquium, ranging from 
metallic point contacts to molecular junctions. 

In this Colloquium we will discuss all these issues 
at the microscopic level. The basic systems we will 
consider consist of a nanoscale junction, namely two 
leads connected by a nanoscale element, with possi- 
bly a third lead controlling some state variable of the 
system, e.g., its local temperature. Typical exam- 
ples are point contacts or quantum do ts placed be- 
tween a two-dimensional electron gas dGodiin et al . 
19991: Ivan Houten et all Il992at iMolenkamp et all 11994 : 



Scheibner et a?. 7 120071 : IStarine et all 1993 ). a molecule 
trapped between a su bstrate and a scanning tunneling 
microscope (STM) tip dBaheti et all 120081: iReddv et all 
2007), metallic wi res (|Ludoph and Ruitenbeekl. 1X9 99) , 
carbon nanotubes (Kim et all l200lb lYu et all I 20051) or 
silico n nanowires (jBoukai et all 120071 : iHochbaum et all 
l2007f) between two metal contacts, etc. Fig. [1] shows a 
schematic representation of the different systems we con- 
sider. The leads are held at different temperatures, which 
allow for the flow of energy (and possibly charge) through 
the junction. Here, we point out that, due to space lim- 
itations, we will not be able to discuss the entire class 
of systems collectively known as "nanomaterials" - com- 
posite layers of various materials fabricated on nanometer 



scales, which show unique electronic properties, often en- 
gineered by adding scattering mechanisms (for instance 
boundary sc attering), that may be beneficia l for energy 
applications (|Maiumdarl . |2004 Is. Vol! 12009) . The inter- 
ested reader may refer to Chen (jChenl . 2005) for systems 
other than those presented here. 

To make the review easier to follow for the reader, we 
have divided it into three main (yet closely related) sub- 
topics. The first one is the transport of heat through 
the system by phonons (lattice vibrations) and electrons, 
which (in linear response) is mainly characterized by 
the thermal conduct ivity k. This is s ue has already been 
revie wed elsewhere ( Galperin et all l2007at IWang et all 
2008), emphasizing the effects of vibrations and focus- 
ing primarily on the method of non-equilibrium Green's 
functions. To make the present review complete, and in 
order to highlight the various theoretical methods and 
the open questions that still pertain to this subject, we 
give it some space here too. In particular, we will dis- 
cuss the different processes that contribute to k and their 
importance in nanoscale junctions. 

The second subject is that of the local temperature and 
heating inside the nanoscale system. This issue is par- 
ticularly subtle, precisely because we are dealing with a 
non-equilibrium process where a temperature difference 
is set at the two sides of the nanojunction. We will ad- 
dress several experimental and theoretical issues and fun- 
damental open questions, such as: How does one define a 
local temperature at the nanoscale in a non-equilibrium 
situation? What determines the local temperature and 
the temperature profile along the system? 

As a corollary of the above studies we are finally led to 
analyze a nearly two-century old and important phys- 
ical law, which so far has eluded a satisfactory theo- 
retical understanding, namely Fourier's law (FL). This 
law, as originally formulated, states that in the pres- 
ence of a temperature difference between the two leads, 
(i) a temperature gradient develops, (ii) the energy cur- 
rent density is proportional to it, and (iii) the constant 
of proportionality is independent of system size. While 
FL was empiricall y postulated fo r bulk systems almost 
two centuries ago (iFourierl . Il822l) and has been derived 
phe nomenological ly for phonons more than eighty years 
ago (|Peierlsl . [l929l ). no simple proof of its validity (or in- 
validity) has ever been derived from first principles, nor 
do we have a well-defined set of conditions to determine 



its validity for a given system (jBonetto et aZ.1 . 120001 ). As 
we will emphasize later, the issue has everything to do 
with the difficulty in defining the basic quantities that 
enter its formulation - namely the local temperature and 
heat current - from a microscopic, quantum mechanical 
point of view. 

The final issue is that of the inter-relation between 
the heat flow and the electron transport through the 
junction, which can be collected under the general name 
of "thermoelectricity" . The central quantity here is the 
thermopower (or Seebeck coefficient) S, which describes 
the voltage drop generated by a temperature difference. 
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A sample of important open questions for this topic are: 
What are the different mechanisms contributing to ther- 
moelectricity? Are they properly taken into account in 
the present theories? What are the state-of-the-art ex- 
periments, and are their results interpreted satisfacto- 
rily? 

All these issues and open questions will accompany us 
for the full length of this Colloquium. We will stress their 
importance for both their fundamental character as well 
as their impact in possible technological applications. We 
will finally point out possible future research directions 
that could explore them in more depth. 

The Colloquium is organized as follows. In Sec. |TT] we 
discuss heat flow in nanoscale systems due to phonons, 
electrons and their mutual interaction, and describe the 
different processes which contribute to it. We review 
both theoretical tools and state-of-the art experiments 
for measuring heat flow in nanostructures. We devote 
Sec. OH] to local temperature effects, and proceed to dis- 
cussing Fourier's law. In Sec. IIVI we discuss thermoelec- 
tric effects in nanoscale junctions. We give a detailed 
account of present theoretical tools, and discuss recent 
experiments, with emphasis on open issues in the field. 
Finally, we conclude in Sec. [V] with some prospects on 
the future of the field. 



II. HEAT CURRENT AND THERMAL CONDUCTIVITY 

Let us start by reviewing the topic of heat cur- 
rent and thermal conductivity of nanoscale junc- 
tions. We will not present full derivations of the 
methods and results. Rather, we will outline only 
the main theoretical tools. The interested reader 
may fi nd extensive a cc ounts i n recent revi ews (Dhar, 
l2008t iGalperin et all l2007bl IWang et all 120081) or 
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book s (jAkkermans and Montambauxl l2007t iDi Ventral . 
2008) where these methods are discussed in detail. In 
addition, we will review recent experimental advances in 
measurements of the thermal conductivity in nanoscale 
systems, with emphasis on the measurement process it- 
self and open questions. 



A. Definitions 

When a nanoscale junction is placed in contact with 
leads held at different temperatures, energy flows through 
it. The original qualitative description for this phe- 
nomenon in b ulk materials is attributed to Fourier 
([Fourier! I1822T) . and amounts to Fourier's law which 
states that a temperature gradient VT induces a ther- 
mal current density linearly proportional to it, namely 



k is the thermal conductivity (such an equation is usually 
valid only in the linear regime). 

In Sees. IIII.BI and IIII.DI we will expand more on the 
significance of the term "temperature" for a system out 
of equilibrium, and its different definitions. Here, we 
anticipate that whenever we do not discuss its mean- 
ing explicitly we call temperature that which is measured 
by a local thermal probe weakly coupled to the system 
and whose temperature has been adjust ed so that t he sus - 
tem dynamics is minimally perturbed ( Di Ventral . [2008) . 
This defin es what we will later call a temperature float- 
ing probe ( Dubi and Di Ventral . l2009d ldl). Note that we 
do not define it in terms of a probe adjusted so that the 
thermal current between the system and probe is zero, 
precisely because we do not have means to measure di- 
rectly the thermal current (although these two definitions 
may give the same quantitative results). In addition, the 
reader needs to keep in mind that while this is an oper- 
ational definition of temperature out of equilibrium, its 
actual experimental determination is far from trivial at 
present. 

The validity of Eq. (JTJ in nanoscale junctions is dis- 
cussed in detail in Sec. IIII.DI Here, we are mainly in- 
terested in the theoretical understanding and measure- 
ment of jth and k, assuming that Fourier's law is indeed 
valid. A relation between the formalism described below 
(Landauer's formula ©) and Fourier's law can be de- 
termined, which requires calculation of thermal conduc- 
tances at larger and larger length scales . Such derivation, 
discussed in other reviews (|Dharl . 120081) implies going be- 
yond the realm of nanoscale junctions and will thus not 
be discussed in detail here. 

It is also convenient to introduce the thermal conduc- 
tance, which is the ratio between the total heat current 
Jth and temperature difference AT = Tr — T^ , 



Jth 

o t h = - hm — — 

AT->0 AT 



(2) 



If the sample is uniform with a constant cross section 
A and length L, the thermal conductance is related to 
thermal conductivity k via o~th = j; K - If t ne sample is not 
uniform, then the relation between thermal conductance 
and conductivity depends on the microscopic details of 
the system. In addition, in analogy with electric circuit 
theory, it is convenient to define the thermal resistance, 
being the reciprocal of the thermal conductance: pth = 

a th- 

Energy can be carried through a nanoscale junction (or 
through a solid) either by lattice vibrations (phonons) 
or by electrons, or both . In insulating bulk materi- 
als the electronic contribution is negligible, while it is 



Jth 



-kVT , 



(1) 



where jth is the heat current density (which may contain 
both phonon and electron contributions, see below) and 



At low temperatures energy can also be carried by the electro- 
magnetic environme nt (photons), an effect which was studied in 
mesoscopic systems {Schmidt et al, 2004) but was not system- 
atically addressed in nanoscale junctions. 
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sizeable in bulk metals. This simple distinction is less 
obvious in nanoscale junctions, where, due to the large 
current densities they can carry 2 , the two contributions 
may be equally important and need to be discussed on 
equal footing. For bulk insulating materials, the theory 
of phonon thermal conductivity based on the B oltzmann 
equation was derive d by Peierls (jPeierlst Il955l ) (see also 
the detailed review (lCarrutherslll96ll) '). The main idea 



is that k is governed by phonon scattering, especially the 
so-called Umklapp scattering (processes that do not con- 
serve crystal momentum), whereby phonons scatter be- 
tween states which are separated (in reciprocal space) by 
a reciprocal lattice vector. 3 Considering a phonon mean- 
free path I (mainly due to scattering by impurities) , sim- 
ple arguments lead to the followi ng relation at high tem- 
perat ures (in three dimensions) ( Ashcroft and Merminl . 
I1976D 



-Ivcv 



(3) 



where v is the velocity of sound and c v is the phonon heat 
capacity at constant volume (in the above equation opti- 
cal and acoustic phonons are considered on equal footing, 
although only the latter ones participate in heat trans- 
port). In a bulk metal a sim ilar relation can be derived 
( Ashcroft and Merminl . Il976l ). where now I stands for the 
electronic mean-free path, c v is the electronic heat capac- 
ity at constant volume, and v is the electron drift veloc- 
ity. Here, a comment is in order. In the case of electrons 
the heat (or thermal) current contains also a contribu- 
tion from the variation of number of particles. In fact, 
let us consider the thermodynamic relation (at constant 
volume) SQ = d£ — /idn, where Q and £ are the heat 
and energy per unit volume, respectively, n is the parti- 
cle number density, and \x the chemical potential. From 
this relation, dividing by the infinitesimal time interval 
dt, we obtain (e is the electron charge) 



Jth = Je 



(4) 



namely, for electrons the heat current has both a contri- 
bution from the energy current, Je, and from the charge 
current J e (there is no such term for phonons, since their 
number is not conserved) . In this review, we will use the 
terms "energy current" and "heat (thermal) current" in- 
terchangeably, but with the understanding that, in the 



2 For instance, in an atomic quantum point contact of a nominal 
cross section of 0.1 ran 2 , to a typical current of 1 /xA corresponds 
a current density of about 10 9 A/cm 2 . This is several orders of 
magnitude larger than in mesoscopic or bulk systems. 

3 For a homogeneous bulk system in which the Umklapp processes 
are suppressed and only "normal" processes occur (namely, pro- 
cesses that conserve crystal momentum) energy can flow undis- 
turbed, giving rise to a diverging k, and such a system cannot 
reach local or global equilibrium. 



case of electrons, one must generally include a contribu- 
tion from the variation of the number of particles (see 
also discussion after Eq. (j6])). 

It is now natural to ask whether these arguments can 
be extended to the regime in which strong material inho- 
mogeneities are the norm, like in nanoscale systems. Be- 
fore we embark in this quest, however, it is worth asking 
why k is such an important quantity in the first place, es- 
pecially since measuring the thermal conductivity at the 
nanoscale is all but trivial. The answer is that k contains 
information regarding two main processes relevant to the 
future applicability of nanoscale systems. The first is the 
rate at which energy is dissipated in and removed from 
the junction. This has an effect on the heating of the sys- 
tem, which may affect its structural stability. The second 
is that k is an important (and limiting) factor in the effi- 
ciency of nanoscale systems as heat- voltage converters (as 
it will be discussed more at length in Sec. lIV[) . Therefore, 
according to the desired use, an ideal nanosystem should 
have opposite thermal properties: for current-carrying 
wires one wishes a high thermal conductance that would 
allow heat to pass through the wire and prevent over- 
heating, and for thermoelectric conversion one requires 
a thermal conductance as small as possible. These re- 
quirements make the understanding, predictability, and 
control of k highly desirable. 



B. Experiment 

In this section we focus on the experimental measure- 
ments of the thermal conductivity in nanoscale systems. 
As already pointed out, a major difficulty in measuring 
k (other than the usual ones related to any measure- 
ments at the nanoscale) stems from the simple fact that 
there is no direct way to measure a heat current. Indeed, 
the only directly measurable quantities are electrical cur- 
rents, voltages and temperatures (the latter also typically 
measured via resistance measurements), and from these 
one deduces n. The main limitation is that the value of 
k as extracted from the experiment may then depend on 
the model used to describe the whole experimental setup 
or device, which may generate some ambiguity. Here, we 
will describe some recent experiments, discuss the meth- 
ods employed in deducing k, and review some of the main 
results. 

A conceptually simple way to measure the thermal con- 
ductance of a suspended nanojunction is the following. 
Consider the schematic system of Fig. [5J The "heater" 
coil is heated by passing a current through it. By measur- 
ing the current and the voltage through the heater coil, 
the power transferred through it is given by the well- 
known relation, P = IV . This power increases the tem- 
perature of the coil to T^. At the same time, the temper- 
ature of the "sensor" coil, T s , is evaluated (by measuring 
its resistance, which is pre-calibrated to correspond to a 
given temperature). If the wire is suspended, then the 
entire heat current should be equal to the power supplied 
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by the heater coil, Q — P, which is related, in linear re- 
sponse, to the temperature difference by 



Q 



-&th(Th — T s ) , 



(5) 



from which the thermal conductance oth can be evalu- 
ated, under the assumption that all the power supplied 
by the electric circuit flows through the junction with- 
out loss, and the thermal conductivity k can then be 
extracted from a microscopic model that relates thermal 
conductivity to thermal conductance (see Sec. Ill.Cj) . If, 
as indeed is the case in many experiment, some of the 
power is lost due to heat diffusion away from the con- 
tacts (e.g. into the substrate) then the Joule heating is 
the sum of the heat flowing away through the contacts 
and that flowing through the wire. 



heater coil 



suspended nano-wire 



sensor coil 




FIG. 2 (Color online) (a) A schematic representation for a 
simple setup to measure the thermal conductance, (b) An 
actual device to measure the thermal condu ctance of boron 
nitrade nanotubes (from (Ch ang et all [2006)). 

This method seems very simple, and was indeed em- 
ployed to measure th e quantum of thermal conductance 
([Schwab et all 1200(3) . However, it needs to be acknowl- 
edged that it has obvious limitations. For one, dissipa- 
tive effects at surfaces or local thermal gradients in the 
heating and cooling parts of the coils 4 may reduce the 
heat flow in the suspended wire. In addition, recent the- 
oretical studies indicate that the contact thermal resis- 
tance between nanowires and substrate plays an impor- 
tant role in determining the overall thermal resistance 
(|Chalopin et all l2008t IZhong and Lukei 120061) . 



More difficult is the determination of the thermal con- 
ductivity k from a model that includes all the effects 
of device geometry and dissipation through the con- 
tacts and substrates. Such models vary for different de- 
vices and geometries ( Chang et aZ.I . [2006. 2008; Shi et all 



I2003D . but share the common feature that thermal con- 
ductances are treated on the same footing as classical 
(charge) conductances, with the same Kirchhoff-like laws 
for the addition of resistances in series (p t h = Y^iPth,U 
with p t h,i the thermal resistance of a single element of the 
circuit) and parallel (p^ 1 = Y^iPthi)- Thus, the mea- 
sured thermal conductivity may depend on the circuit 
model used, which makes it hard to compare between 
different experiments. This means that when perform- 
ing a measurement, one is in fact measuring the ther- 
mal conductance of the system of interest embedded in 
that specific device. Nevertheless, this method was used 



structures, mainly carbon nanotubes 


Brown et all 2005 


Chang et al, 20061 20081: Chiu et al. 


. 2005; Fuiii et al. 


2005; Kim et al. 20011 Shi et all 2003: Yu et all 2005) 



but also nanotubes of other material s 
120081: IChen et all 120081: iLi et all I2003D . 



Chang et al. 
Some experi- 



mental fe atures are universal , like ballist i c ther mal con- 
ductance ( Brown et all [20051 : IChiu et al. , 2005), a value 
of k which is orders of magnitude larger than the 
bulk value for carbon nanotubes (~ 3000W/K at room 
temperature), an increase of thermal conductance with 
nanowire diameter, or a peak of the thermal c onduc - 
tance at - 320K ([Fuiii et all 120051: iKim et all 1200 if) , 
attributed to the onset of Umklapp phonon scattering 
processes. However, other features, such as the detailed 
power-law dependence of n on temperature vary between 
experiments, indicating that this is not a universal fea- 
ture, and depends on the details of the experimental 
setup. 

Other experimental approaches to measure n have 
bee n introduced in t he literature. For instance, Pop et 
al. (|Pop et all 120061 ) have used high currents to induce 
heating in a single- walled carbon nanotube, with a model 
to relate the current-voltage (I-V) characteristics to the 
high-temperature thermal condu ctance. In another ex - 
ampl e, the so-called 3w method (|Cahilll. 19901 Lu et al. 
2001), was used to study na notubes ([Bourgeois et al. 



4 Recall that one can destroy and create phonons at the surfaces 
of a material. 



20071: IChoi et all 120051 12006D . In this method, an a.c. 
current is applied to the sample which also acts as a 
heater. From a simple derivation one finds that the third 
harmonic of the voltage drop across the sample is re- 
lated to the thermal conductivity of the sample (at small 
frequencies of the current). Using this method, the au- 
thors found a deviation of the thermal conductance from 
a cubic dependence on temperature for Si nanowires, in- 
dicating a dimensional crossover at low temperatures. 
Both these methods rely on current-induced self-heating 
of the sample (rather than direct heating by an external 
source). In a third example, laser- induced heating and 
Raman spectroscopy (already use d in various nanoscale 
systems such as graphene ribbons ( Balandin et all 120081: 
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ICalizo et al. I, Hqq3)) have been used to determine the 
local temperatures (|Deshpande et all l2009t iHsu et all , 
2009) and extract the thermal conductance of carbon 
nanotube bundles. The main disadvantage of this 
method is that to obtain the thermal conductance one 
needs to assume a value for the optical absorbtion of the 
sample, which is usually unknown. 



C. Theoretical Methods 

We now provide a brief description of the theoretical 
methods most commonly employed to describe energy 
flow, with an eye on their strengths and limitations. 



1. Single-particle scattering approach 

Many theoretical calculations of thermal conduc- 
tance a re based o n an a pproa ch pioneered by Lan- 
dauer (|Landauerl . I1957L Il970h in the context of 
charge trans port in mesoscopic and nan oscop ic sys- 
tems dDattal . Il997t iDi Ventral . 120081: llmrvl Il997h . The 
same ideas have been generaliz ed to phonon transport 



where 



through a nanoscale junction (lAngelescu et al\, 
Blencowel fl999l 12004 (l5har and Rovl. 120061: iRegol 



1998; 



2001 



Rego and Kirczenow . 19981 : Segal et al. . 2003 ). 

The basic tenet of this approach is that one assumes 
the leads non-interacting (ot herwise no closed form for 
the current can be obtained (|Di Ventr so that 

a convenient basis, such as plane- waves, can be chosen to 
develop state vectors for both types of particles, either 
phonons or electrons. As a further conceptual simplifi- 
cation, the leads are thought to be adiabatically "con- 
nected" to reservoirs whose only role is to define the 
occupation of the scattering states according to a local 
equilibrium Bose-Einstein (BE) distribution for phonons 
or a Fermi-Dirac (FD) distribution for electrons. Once 
this occupation is set, the particles are free to propagate 
in the leads before scattering at the lead-system inter- 
face. Charge and/or energy current is then determined 
by an electrochemical potential difference and/or a tem- 
perature difference between the reservoirs. 

Most of the calculations also assume that the particles 
in the sample are either truly non-interacting or inter- 
acting at a mean-field level (which is the same from a 
formal point of view) . In this case the current is simply 
proportional to the probability for the particles to cross 
the sample from one electrode to the other. For instance, 
in the case of phonon transport, phonon states at a given 
energy huj, scatter off the junction and may be cither 
transmitted through it or reflected back. The probability 
to be transmitted through the junction is characterized 
by the transmission coefficient T(w). The expression for 
the heat current is then simply 



f°° du; w 
= J 2^ nijJ T{Lu){g L - g R 



(6) 



9l{r) 



L/ (exp ( 



k B T 



L(R) 



1 ] are the distribu- 



tion functions of phonons in the left (right) lead. From 
Jth one can then evaluate the thermal conductance ac- 
cording to Eq. ([2]). 

Within this approach the electronic contribution to the 
heat current is calculated similarly, where in Eq. ^ one 
makes two changes, namely (i) the BE distribution func- 
tions are replaced by FD distributions, and (ii) the en- 
ergy in each reservoir is measured from the respective 
electrochemical potential, and /i^ for the left and 
right reservoir, respectively, i.e., hu> — > hu> — [1l,r (see 
Eq. (|4])). In linear response this leads to the substitu- 
tion hu — > huj — (fiL + A 1 i?J/2 in the energy term that 
multiplies T(uj) in Eq. ©. 

To actually evaluate ath, one has to compute the trans- 
mission coefficient T(w). To this aim several methods 
hav e been employed , such as the use of continuum mod- 
els (lAngelescu et q/1 Il998l) . boundary condition method 
dWang and Wangt 1 2006) , mode - matching method (lAndol . 
119911: iKhomvakov and Brocksl . 12004 iTing et all Il99a) 
and scattering or t ransfe r matrices (|Di Ventra and Land . 
120021: iTong et all Il999h . All these methods are fun- 
damentally equivalent, and in fact have their origin in 
the single- particle elastic sc attering theory of conduction 
(see, e.g., (|Di Ventral . 1200 81). whereby one can write the 
transmission coefficient as a sum of all the partial prob- 
abilities of transmission Tif(ui) from one of the momen- 
tum states of the incoming (i) particle (whether elec- 
tron or phonon) at energy Huj to one of the momentum 
states o f the outgoing (f) pa rticle at the same energy, 
namely ( Buttiker et aLl . ll985l) 



(7) 



where r is a sub-matrix of the scattering matrix with 
dimensions Nr x Nl, with Nr and Nl the number of 
channels in the right and left leads, respectively, at en- 
ergy hjj. This result can be cast in another equiva- 
lent form in terms of single -particle Green's functions via 
(|Meir and Wingreenl fl992l) 



T{u) = Tr{G r T L G a T R } . 



(8) 



where G r ^ is the retarded (advanced) single-particle 
Green's function corresponding to the interaction of 
a "central" part of the junction with the elec- 
trodes and ^l(r) describe the "rate" at which par- 
ticles scatter between the leads and the central 
part of the junction. It has been r e-deriv e d for 
thermal trans p ort by many auth o rs dDhar . 2008 : 



2003; Ozpin eci and Cirac 
Wang et all l2006l |20(T 
200 



Uracil. 1200: 
20071: lYanl 



Galperin et all l2007al iMingol 120061: iMingo and Yang 



200 11: ISegal et all 12 003: 

Yamamoto and Watanabe 



Arguably the most universal result obtained from the 
Landauer formula ([5]) is that of thermal conductance 
quantization. Similarly to the quantization of electrical 
conductance in ideally one-dimensional (ID) electronic 
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systems ( van Wees et q/.l . ["l988l) . at low temperatures the 
thermal conductance (per phonon mode) was predicted 
to acquire a quantized value 



ctq 



n 2 k 2 B T 
3h 



(9) 



where h is Planck's constant ( Greiner et all 



Maynard and Akkcrmans, 
Rcg o and Kirczenowl Il99 



119851: 
This 



Pendrv 



1997 



1983; 



result is readily 
derived from Eq. [5] in linear response by setting the 
number of modes to unity, and letting the transmission 
coefficient to be one, i.e., T(u>) = 1. 

The fact that this conductance is material-independent 
relies on the fact that, like in the electronic case, in ID 
the phonon density of states is exactly proportional to the 
inverse of the group velocity. Remarkably, thermal con- 
ductance quant ization does not depend on the statistics 
of the carriers (|Regol . l200lD . Indeed, op wa s expe rimen- 
tally measured for phonon s (ISchwab et all l2000l ) , elec- 
trons (jChiatti et all 120061: iNicholls and Chiattil . l2008h 
and even photons (|Meschke et all |2006[)~ 

Another application of the Landauer formula © 
has been in the study of geometrical and tempera- 
ture effects on thermal transport. To give a few ex- 
amples, this approach has been used to understand 
the role o f defects on the t hermal conductance of a 
nanowire (|Chen et all l2005aft . the effects of different 
geometries such as stubs, T-j unctions and concavities 
(|Peng et all l2007t iTang et all 120061: IXie et all l2008j) . 



periodic m odulations (|Tang et all 120071). and surface 



rough ness (jKambili et all Il999l : ISantamore and Crossl 
1200 ll ). As a general rule, disorder and temperature are 
found to have competing roles: disorder tends to reduce 
the thermal conductance (by decreasing the transmission 
coefficients of the different transport modes), and a tem- 
perature increase usually results in a larger thermal con- 
ductance, due to an increased number of modes which 
participate in the thermal transport. 

The interplay between the two processes can result 
in i nteresting phenomena . For instance, Santamore et 
al. ( Santamore and Crossl . l200lh showed that disorder 
in the form of surface roughness may generate a non- 
monotonicity in crth with increasing temperature, with a 
slight decrease (below the quantum of thermal conduc- 
tance) followed by a rise of oth with increasing tempera- 
ture, in similarity to t he exp erimental results of Schwab 
et al. (jSchwab et all [2000) . Their results (shown in 
Fig. [3]) are explained as follows: at very low tempera- 
tures, there is only one mode which contributes to the 
thermal conductance. As temperature increases, scat- 
tering of that mode off the surface roughness increases, 
generating a decrease in the thermal conductance. As the 
temperature is raised even higher, higher modes start to 
participate in the thermal transport, giving rise to an 
increase in the thermal conductance. 

This ties with the use of the scattering approach 
to thermal conduction in real materials, which comes 
about from using realistic phonon spectra (e.g., as ob- 




0.4 ftB 
Temperalure T (K) 



FIG. 3 Thermal conductance of a quasi-lD wire with sur- 
face roughness, exhibiting inhomogeneous thermal conduc- 
ta nce. Points correspon d to experimental data of Schwab et 
al. dSchwab et all [20 00 ) , and the solid line is the theoretical 
curve. From (Sant amore and Crosd I2001I ). 



tained from experiment or first-principles approaches) 
in combination with ground-state density-functional the- 
ory (DFT) calculations to obtain the scattering coef- 
ficient T(uj). To give several examples, Tanaka et 
al. have combined geometrical structure (i.e., realistic 
shape of the wire) with real material parameters to 
study the onset of the thermal conduct ance quantiza- 
tion in GaAs and silicon nitride wires ( Tanaka et all 
2005). The thermal conductance of nanowires made 
of, e.g., Si, Ge and Ga A s wa s studied by s everal 
authors (IMingo arid Yangl 120031: iMingo et all 120031: 
iTanaka et all \200A IWang and Wangt 120071) . Much at- 
tention has been given to carbon-based struc tures, such 
as ca r bon nanotubes , grap h ene and graphite (lLan et al. . 
20091: iLii and Wand. 120081: IMingo and Broidol. l2005aUb 



Yamamoto and Wat anab e 



2006: Yamamoto et aZ.I .[2004: 



Zhan g and Li . 20051 : Zimmermann et al. . 2008 ). Another 



exampl e is the recent study of isotope and disorder 
effects ( Murphy and Moore] . 120071 ) , specifically in ca r- 
bon and boron-nitri de nanotubes ( Savic et~al\ . l2008aU bl; 
IStewart et~al\ . l2009h . 

Some universal conclusions arise from these calcula- 
tions. For instance, a dimensional crossover from three- 
to one-dimensional transport (manifested by, e.g. a 
change in the temperature dependence of the thermal 
conductance) occurs in many systems as the diame- 
ter of the nanotube decreases, the length scales deter- 
mined by the wavelengt hs of the typical phonon modes 
(|Wang and Wand . l2007D . Also, disorder in various forms 
(local defects, surface roughness, etc.) has a dramatic 
effect on the thermal conductanc e, as it influences th e 
scattering of the different modes ( Rov and Dharl . l2008h . 
Due to the translational invariance of the lattice, long 
wavelength (or zero-frequency) modes are always con- 
ducting, while short wavelength modes are scattered by 
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disorder. Since the short wavelength modes participate 
in the thermal transport only at high temperatures, it is 
found that the low-temperature thermal conductance is 
less affected by disorder and defects. Finally, the thermal 
conductance of molecul ar junct ions has also been widely 
studi ed (jGalnerin et alll2007allMiiigol . l200fit lSe gal et all . 
2003). It is found to be strongly dependent on a multi- 
tude of factors, among which the phonon spectrum of 
the molecules, the degree of localization of the molecular 
modes, the molecule-lead coupling, non-harmonicity (i.e. 
phonon interactions), etc. 

It is important to stress once more that Eqs. ©, Q), 
([8]) , and indeed the whole Landauer approach, are based 
on some strong assumptions, which may breakdown in 
nanoscale junctions and under certain experimental con- 
ditions. The first assumption is that the system is 
"closed" , in the sense that it does not dynamically inter- 
act with its environment. The latter only provides the 
boundary conditions and the relevant parameters (like 
the temperatures, chemical potentials, etc.). The second 
assumption is that the leads are ideal, i.e., are unaffected 
by the proximity to the junction (either in their struc- 
ture or in the distribution of particles) and support well- 
defined single-particle states. In addition, it is assumed 
that "dissipation" takes place at the (infinitely far) edges 
of the leads and that the temperature (and chemical po- 
tential for electrons) is uniform in them. Most critically, 
the approach does not provide information on the dy- 
namics of the system. Therefore, transient, memory and 
non-linear dynamical phenomena are beyond its reach. 
A further issue arises when one uses ground-state DFT 
in combination with the Landauer approach: one is effec- 
tively using a ground-state theory for a non-equilibrium 
problem. This issue cannot be solved by knowledge of the 
exact ground-state exchange-correlation functional, and 
as such, the use of ground-state DFT in this context can 
only be viewed as a mean-field a pproximation. This has 
been explicitly demonstrated in (jVignale and Di Ventral . 
2009), where for the case of electrical conductance it was 

•H- 

shown that the exact resistivity tensor p can be written 
as 



P = Ps + Po 



(10) 



where P s is is the resistivity tensor of a noninteract- 
ing system in the presence of a static potential V s 
that reproduces the exact ground-state density, and 

p xc is a dynamical contribution related to dynamical 
exchange-correlation effects, and which does not van- 
ish even in the zero-frequency (d.c.) limit. A possible 
way out would be to use a fully dynamical approach 
(e.g., the microcanonical picture of tr ansport as sug- 
gested in (|Di Ventra and TodorovLl2004)) c ombin ed with 
time-dependent DFT ( jRunge and Gross! Il984l) . This 
appro ach (recently imple mented to study charge trans- 
port (|Cheng et all l2006h ) would provide, in principle, 
the exact thermal total current, if the exact dynamical 
exchange-correlation potential is known. However, we 



are not aware of any calculation of thermal current along 
these lines. 



2. The role of interactions 

Up to this point the system Hamiltonian has been 
assumed to describe single particles with interactions 
included at most at the mean-field level. As briefly 
mentioned above, many-body correlations can be ac- 
counted for within a time-dependent DFT approach, 
namely within an effective single-particle picture. Al- 
ternatively, the effect of interactions beyond mean- 
field, could be explicitly included via the so-called non- 
equilibrium Green's functions formalism (NE GF) (see, 
e.g., N. Mingo, chapter 3 in (|S. Volzt 120091) 1 In this 
approach one solves equations of motion for appropriate 
single-particle Green's functions th at can be conveniently 
defin e d on the Keldysh contour ( Kadanoff and Bavml . 
Il962t iKeldvshl . 11964 ). In its exact formulation, the 
NEGF has however limited practical utility, since if one 
assumes particles interacting - beyond mean-field - in 
the whole system (leads plus nanostructure) no closed 
equation of motion for the single-partic le Green's func- 
tions can be obtained l|Di Ventral . 120081 ). Instead, it is 
common to assume (as in the Landauer approach) that 
the leads contain non-interacting particles and interac- 
tions are confined within a "central" region containing 
the nanostructure. This is a strong assumption and may 
not always correspond to the physical problem at hand 
and/or its experimental realization. 

If one makes the assumption of non-interacting parti- 
cles in the leads, and assumes that a steady-state has 
been reached in the long-time limit (not an obvious 
statement either), the equation of motion for the differ- 
ent single-particle Green's functions can be closed and 
the NEGF provides a compact expression for the to- 
tal current similar to that derived for electron transport 
( Meir and Wingreenl Il992l) . given by 



1 r°° 

J = — J hLodu[{G r -G Q )(£<-£<)+ 

+^G < (^ ^^ -^ L )] , 



(ii) 



where G r,a,< are the retarded, advanced and "lesser" 
single-particle Green's functions, respectively; are 
the "lesser" self energies of the a = L,R leads and 
r a = — £°) namely the difference between "re- 

tarded" and "advanced" self energies (the explicit in- 
dependence of all these quantities has been omitted). 
The first term on the right hand side of Eq. (jlip may 
be interpreted as describing the current from the bias- 
induced difference in the coupling to the leads, while 
the second is related to the non-equilibrium distribution 
function in the interacting region. The single-particle 
Green's functions can represent either phonons or elec- 
trons, and should be calculated in the presence of interac- 
tions. In the mean-field approximation, Eq. (|11[) reduces 
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to Eq. ([6]) (or its equivalent form for fermions). Many- 
body perturbation expansions to compute these Green's 
function s have been performed for simple model Ha mil- 
tonians ( Galperin et all [2007at Eu" and Wand . 120071) but 
it is no easy task to introduce interactions (beyond mean 
field) in realistic systems. 

The NEGF could also be used to study the effects of 
electron-phonon interactions. In that case as well, how- 
ever, quite strong approximations need to be made in or- 
der to have an analytically tractable theory. For instance, 
if one assumes electrons interacting with each other at 
a mean-field level, but interacting in a "central" region 
with non-interacting phonons, the heat current can be 
approximated as a sum of contributions from both elec- 
trons and phonons, J = J e i + J p h , each component cal- 
culated with the help of Eq. (TTTI) . The key ingredient 
here is that, due to the electron-phonon interaction, the 
self-energy of phonons includes an electronic contribution 
and vice versa. These contributions can be calculated in 
a perturbative way. However, this is clearly an ideal- 
ization, since it neglects correlated electron-ion motion, 
which, in principle, does not even allow the total thermal 
current J to be separated into two distinct contributions 
from the two particle species. Along the same lines of 
reasoning, the effects of phonon-phonon interaction have 
been studied (|Liu and Yil . 120061: Wined 120061: IXu et all 
|2008[ ). According to these results both electron-phonon 
and phonon-phonon interactions decrease the thermal 
conductivity. However, we need to stress once more that 
due to the large current densities nanoscale systems carry 
- and hence the large number of scattering events per 
unit time and unit volume - it is not a simple task to in- 
clude all the relevant physical scattering mechanisms in 
the present non-equilibrium case. An example of this is 
the possibility of phonon modes in the junction which are 
weakly coupled to the bulk modes of the electrodes. In 
this case, these "localized" modes may be energetically 
"pumped" by scattering with electrons or other phonons 
before energy could efficiently be dissipated away. This 
physical situation is beyond second-order perturbation 
theory and more work in this direction is thus highly 
desirable. 



3. Molecular dynamics 

Another method to evaluate the thermal conductivity 
which is gaining increasing popularity is that of molecu- 
lar dynamics (MD). Basically, molecular dynamics comes 
down to solving the classical equations of motion of the 
system numerically. The origin of the method in the 
present context can be tr aced back to the se minal work 
of Fermi, Pasta and Ulam (jFermi et oi.l . ll955l ). where the 
energy transfer in non-harmonic lattices has been stud- 
ied numerically. Since then it has been wi dely used to 
study heat transp ort in classical ID systems ( Dharl . [2008t 
ILepri et~al\ , 120031) . It has also been generalized to study 
quantum effects, by providing appropriate boundary con- 



ditions ( Wang et l2008h . These approximations, how- 
ever, should be thought of as quasi-classical, since the 
microscopic dynamics of the system is described by clas- 
sical Newtonian equations of motion, and the quantum 
nature is only introduced via indirect conditions (such as 
the noise in a Langevin term) . A big advantage of molec- 
ular dynamics is the ability to model realistic systems and 
geometries in a rather straightforward way. The forces 
between atoms are evaluated from realistic parameters, 
so that different geometries, impurities, structures, etc. 
are easily taken into account. 

In order to calculate the heat transport directly from 
MD, one needs to account for a finite temperature in 
the system. This is usually done in linear response 
by adding to the Newtonian equations of motion a 
Langevin fluctuating term which satisfies the fluctuation- 
dissipation relation, i.e., the two-time correlation func- 
tion of th e current is proport ional to the temperature 
(see, e.g., (|Van Kampenl . 120011) ). Alternatively, a Nose- 
Hoover thermostat is introduced, in which a fictitious 
coordinate is adde d to the r eal coordinate to maintain a 
finite temperature (|Hooverl . Il985t iNosel . 1 19841) . 

Once a finite temperature is set, there are two main 
methods to calculate the thermal conductivity. The first 
(sometimes called equilibrium MD) is via the linear- 
response Green-Kubq formula (|Dharl. 120081: ILepri et all 
120031 : iLuttineen . I1964D 



Wk B T 2 



(Jth{t)J th (P))dt 



(12) 



where V is the volume, k B the Boltzmann constant, T 
is the system temperature, Jth(t) = J drj t h(r,t) is the 
integral of the heat current density, j t h(r,t), over the 
entire system, and the brackets denote equilibrium en- 
semble averaging in the absence of a thermal gradient. 
However, the Green-Kubo equation has two main weak- 
nesses. The first is that it is derived in the thermody- 
namic limit a nd therefore its use in finite systems is not 
well justified (|Kundu et al .1. 120091) . Secondly, one needs 
to assume that a small temperature gradient (the exter- 
nal perturbation) ensues in the system, which may not 
be the case in every experiment. However, its relative 
simplicity makes it a good starting point in many cases. 

An alternative method (also known as nonequilibrium 
MD), still based on molecular dynamics, is the one in 
which the system is held in contact between two heat 
baths of different temperatures. Once the dynamics 
reaches a steady state, the temperature profile and the 
local heat currents can be calculated, from which the 
thermal conductivity is extracted. Here lies one of the 
disadvantages of the model, since the definition of the 
local heat current requires defining a local e nergy opera- 
tor, which is not always a unique quantity ( Lepri et all 
120031 : IWu and Segal 120091 ). Likewise, a local temper- 
ature needs to be defined and evaluated; a somewhat 
tricky issue to which we will come back in Sec. IIII1 At 
high temperatures (where the distribution function is 
practically classical and quantum effects are negligible; 
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say at temperatures higher than the typical vibrational 
mode temperature) one may define the local tempera- 
ture as the kinetic energy of the atoms (via the equipar- 
tition function), but this assumption breaks down at low 
temperatures, and one needs to use a definition of tem- 
perature which rests on th e equilibrium distribution of 
phonons (|Wang et al. 1. 12008b . This yields a quasi-classical 
treatment (which is somewhat better than a fully clas- 
sical treatment at low temperatures), but leans on the 
assumption that the phonon distribution resembles its 
equilibrium form, which may not be the case in this non- 
equilibrium problem. On the other hand, the obvious 
advantage of this method is that it does not rely on 
any thermodynamic-limit assumptions and is thus ap- 
plicable for any system size, which is important for the 
study of realistic nanosca le systems. For instance, Yang 
et al. ( Yang et all l2010h recently used the method to 
study Fourier's law and thermal conductance of realistic 
Si nanowires, and showed that Fourier's law breaks down 
in these systems (see Sec. IIII.Pp . Studies along simi- 
lar lines have been recently performed to i nvestigate the 
thermal conduc t ance o f carbon nanotubes dBerber et all 
\200dt iHu et all 120081: iPadgett and Brennerl. |2004. Si 



wires (Henry and Chenl . I2008bt iPonomareva et all , 2007 



Yang et al. . 2008h . diamond nano-rods (Padgett et al 



20061 ) andpolvethvlene chains (|Henrv and Chen . 2008a), 



to name only a few recent studies. 

An additional method, related to MD, is that of lattice 
dynamics models. In this method the phonon dispersion 
relations are obtained by calculating the direct change 
in energy due to atom displaceme nts, using force fields 
obtained from DFT calculations flFeldman et al. |, l200d: 
iRen et all l2006t iTurnev et all . l2009h . 

The abundance of literature makes it hard to de- 
scribe universal features of the thermal conductance, 
which seems to strongly depend on the details of the 
model and/or material. Specifically, a t h is very sensi- 
tive to the phonon s pectra and to phonon loca l izatio n 
(|Dhar and Lebowitj . l200l IZhernov and Chulkinl . l2000h . 
which are in turn sensitive to material, geometry and dis- 
order, surface roughness, and more. The rationale behind 
these studies is that by uncovering the detailed influence 
of these parameters on a t h , theory may provide guidance 
to experiments and even suggest new materials with op- 
timized thermal properties. 



III. LOCAL TEMPERATURE AND HEATING 

A. General remarks 

When a current passes through a classical resistor, the 
latter heats up. This phenomenon is known as "Joule 
heating" . It is a consequence of the inelastic relaxation 
of electrons in the r esistor which transfer e nergy to the 
surrounding lattice ( Ashcroft and Merminl . Il976l ). In a 
nanoscale system, such as a molecular junction or an 
atomic wire, electrons can analogously scatter inelasti- 



cally off the phonons (i.e., the vibrational modes of the 
structure). However, since electrons typically spend very 
little time in the junction region, one might naively think 
that their inelastic scattering rate is negligible with con- 
sequent little heating of the junction itself. This con- 
clusion, which is, for instance, at the heart of the Lan- 
dauer scattering approach where all dissipation is as- 
sumed to occur only in the "reservoirs" , does not take 
into account the fact that due to the small cross-section 
of nanoscale systems, the current density at the junc- 
tion can be very large (typically much larger than in 
mesoscopic and macroscopic systems, see footnote [2). 
This implies that the power per atom in the junction 
can be very large, po s sibly leading to large local heat- 
ing (IChen et all I2003L l2005ri iTodorovl Il998). The rate 



at which this power is then dissipated back to the elec- 
trodes determines the effective local (and out of equilib- 
rium) temperature of the junction. 

In addition, current-carrying electrons can transfer 
energy, via inelastic electron-el ectron interactions, t o 



other electrons in the Fermi gas ( D'Agosta et all l2006h . 
This effect is generally small in macroscopic systems. 
However, similarly to the increased rate of electron- 
phonon scattering in nanoscale junctions due to the 
large current densities, the inelastic scattering rate of 
electron-electron interactions may increase in nanoscale 
syst ems leading to a local heating of the electron liq- 
uid ( D'Agosta et all l2006f ). This effective higher tem- 
perature of the electrons may influence the local ionic 
heating due to electron-phonon interactions and thus can 
be indirectly measured, by measuring the local tempera- 
ture of th e ions or the broadening o f inelast ic conductance 
features (|D'Agosta and Pi Ventral l2008a|) . 



An obvious reason why local temperatures and heating 
are such important phenomena lies in the fact that sub- 
stantial heating of a nanoscale system leads to the system 
instability and eventually to the breaking of atomic bonds 
(ITeramae et al. L l2008tlTsutsuT et al. l. l2008bHWard et all 
2008]). A different and even more fundamental interest in 
these phenomena arises in the context of Fourier's law, 
Eq. ([1]), that we will discuss in Sec. IIII.PI Of course, 
at the nanoscale, it seems inappropriate to discuss the 
scaling of the thermal conductance with length, since 
this is an as y mpto tic (in terms of system size) property 
( Lepri et~al\ . [2003). Thus, one is left with the simple 
question: under which physical conditions does a uni- 
form temperature gradient develops in a nanoscale sys- 
tem held in contact between two heat baths of different 
temperatures? 

In this section we discuss all these issues. We re- 
view the various mechanisms which give rise to heat- 
ing in current-carrying junctions, using simple arguments 
and models, followed by some basic results obtained from 
more elaborate models. We then turn to discuss the on- 
set of a temperature gradient, analyzing a molecular wire 
junction in terms of the theory of open quantum systems, 
discussed in some detail in Sec. IIV.C.21 
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B. Heating in current-carrying nanostructures: theory 

1. Various definitions of out-of-equilibrium temperature 

In order to discuss local heating, the first question 
one should ask is: how is a local temperature defined 
and calculated? Since temperat ure is a thermodynamic 
quan t ity, some caution is need ed ( Hartmann and Mahlerl . 
120051 : iHartmann et all , l2004afbh . Apart from the de fini- 
tion of temperature that we have given in Sec. III. At and 
which we will use also in Sec. IIII.Di we here report sev- 
eral other notions of local temperature (not necessarily 
leading to the same quantitative results) and their micro- 
scopic origin, which were used to study local ionic heating 
in atomic junctions, each with its own pros and cons. 

Kinetic definition - An intuitive definition of local tem- 
perature is to relate it to the local kinetic energy of 
the ions, i.e. (imt) 2 ) ~ 3fceT/2. However, this defi- 
nition, mainly used in molecular dynamics simulations 
(see Sec. III.C3[) . has several drawbacks: (i) it relies on 
the equipartition theorem which is strictly proven in the 
thermodynamic limit only for systems whose energy is 
quadratic in the particle momenta (as for non-interacting 
systems), and does not encompass any quantum effects, 
(ii) One needs to define an average kinetic energy over 
some length scale, while the quantum nature of particles 
may preclude such definition. 

Local phonon mode - Consider a phonon mode some- 
how coupled to the system and vary its temperature in 
such a way that no heat flows between that mode and 
the system. This idea is somewhat similar to the idea 
of connecting an external bath to a system and impos- 
ing that no heat current flows between the system and 
bath, which was suggested to study the onset of Fourier's 
law i n one-dimensional systems, b o th classical and quan- 



Bonetto et all . l2004t iDharl . 120081 : iDhar and Rovl . 
Rovl . 20081) . This idea was recently used to study 



turn 
120061 : 

the local temperature of a model molecu l ar junc tion us- 
ing the NEGF formalism (|Galperin et aUl2007allbl ). The 
main result is the existence of two voltage thresholds. 
The first is at the voltage which corresponds to the vi- 
brational energy of the phonon, eV ~ hujQ, at which lo- 
cal heating starts to occur and the temperature increases 
abruptly. The local temperature then remains roughly 
constant, until it rises again when the bias is so large as to 
encompass the molecular conduction window (i.e., both 
the HOMO and LUMO states). The disadvantage of this 
method is that the temperature of the mode depends on 
the microscopic details, i.e., the phonon excitation energy 
hujo and/or the electron-phonon coupling. 

Distribution function definition - A slightly different 
model of local temperature is to connect a phonon mode 
to the nanoscale system, but instead of determining its 
temperature self-consistently, its distribution function is 
compared to an equilibrium distribution function with 
a given temperature, which is tuned to give the best 
comparison. Clearly, the disadvantage of this method 
is that the non-equilibrium distribution function may 



be very different from th e equilibrium one ( Koch et all 
120061 : iPekola et all , |2004 ). The last two methods were 
compared, and were found to give similar local temper- 
atures at large bias (compared to the typical vibrational 
modes, implying strong non-equilibrium and population 
of higher modes), but deviated from each other substan- 
tially at low biases. In fact, the second method turned 
out to give erroneous results in the zero-bias limit, when 
one expects the temperature to be the same as that of the 
leads. This is precisely because an equilibrium form for 
the phonon mode was assumed, although even with no 
current the distribution function of the mode may have 
contributions arising from the coupling to the electronic 
(and other phononic) d egrees of freedom in the junction 
|Galperin et aUl2007aft . 

Definition from dissipated power - A microscopic the- 
ory which relies on first-principles was sug gested by Chen 
and coauthors (jChen et all 120031 l2005bf ). The method 



is as follows. As a starting point, the electronic scatter- 
ing states are calculated using ground-state DFT. The 
electron-phonon coupling for the different modes is also 
calculated using first-principles approaches. Using per- 
turbation theory then one can calculate the power dis- 
sipated into the junction from current-carrying states. 
This power is then compared to the rate at which heat 
escapes the junction, typically assumed as ct^T 4 with 
dth the thermal coefficient that can be estimated from 
a microscopic model, and T the effective temperature 
of the junction (jChen et all 120031 l2005fj) . A result of 
these calculations is presented in Fig. [U where the local 
temperature as a function of bias was calculated for a 
benzene-dithiol (BDT) junction and a gold-atom point 
contact. The results indicate that, for a given bias, the 
BDT junction heats up less than the gold-atom junc- 
tion, due to better thermal coupling with the electrodes 
and larger resistance to electrical currents (see Eq. ITS!) . 
This result is also confirmed by experiments on sim- 
ilar systems ( Huang et al] , l2006t iTeramae et al. 



liar systems (Irluang et all UUUbt lleramae et all I^UUet 
iTsutsui et all 120071) . While not visible from Fig. gj 



theoretical results of the threshold voltage for heating 
- see Eq. [13] - are also in good a greement with exper- 
iments l)chen et all 120031 l2005fj) . The same method 
was used to stu dy local heating in alkane chains of dif- 
ferent lengths (jChen et al I l2005tJl . It was predicted 
that, at fixed voltage, heating decreases with increas- 
ing chain length, which is due to increased resistance 
to el ectron flow; a resu lt also confirmed experimen- 
tally ( Huang et aZ.I . l200"7h . More recently, the same ap- 
proach was used to study the effect of different isotope 
substi tutions on t he heating in hydrogen molecular junc- 
tions (IChenl . 120081) . It was found that local heating is very 
sensitive to isotope effects since the electron-phonon cou- 
pling constant is inversely proportional to the ionic mass. 

The method described above has the advantage that 
it can treat realistic systems. However, its main draw- 
back is that it relies on the assumptions of the Landauer 
approach - see Sec. IIV.C.1I - and its practical implemen- 
tation employs ground-state DFT, which, as we have dis- 
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FIG. 4 (Color online) Local temperature as a function of bias, 
calculated from a scattering theory approach, for a benzene- 
dithiol molecular junction ( dashed line) and a gold-atom point 
contact (solid line). From (|Chen et all\2063 \. 
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FIG. 5 (Color online) A schematic representation of the 
mechanism of ionic heating in nanoscale junctions. The elec- 
tric current dissipates a fraction aV 2 /R of its power in the 
junction, depending on the strength of the electron-phonon in- 
teraction. This power is then dissipated to the phonons in the 
electrodes in the form of a heat current. The balance between 
the power flowing into the junction and the heat current Iq 
flowing out of the junction determines the ionic temperature 
of the junction. 



cussed at length in this review, does not take into account 
properly all dynamical effects. 



2. Ionic heating 

After discussing various definitions of local tempera- 
ture, we are now in a position to discuss local heating. As 
described previously, we consider here a junction, com- 
posed of leads (which are assumed to be held at local 
equilibrium), and a nanoscale system which has both 
electronic and vibrational degrees of freedom. Even in 
the presence of current, we can assume that in the leads, 
far away from the nanojunction, electrons and phonons 
reach the same temperature Tq . In the junction, how- 
ever, the electrons and phonons may have different tem- 
peratures, T e and T p h, respectively. These temperatures 
depend on bias, strength of electron-phonon and electron- 
electron interactions, the coupling of phonons with the 
bulk phonons in the leads, as well as the transmission 
properties of the electrons. 

Let us start by discussing the temperature of the ions 
in the junction, or the phenomenon of local ionic heat- 
ing (see schematic in Fig. [5]). We start from some sim- 
ple considerations assuming first no i n elastic electron- 
electr on scattering occurs (|Di Ventr 1 120081: iTodorovl 
1998). The power of the entire circuit (nanojunction 



5 The extent to which this statement is correct depends on the 
current density in the leads. If this current density can be as- 
sumed to be zero, then the leads are at an ideal global thermal 
equilibrium, with electrons and ions sharing the same temper- 
ature. Otherwise some difference (albeit extremely small) may 
arise between the lead temperature of the ions and electrons. 



plus power source) is given by V 2 / R, where V is the 
source bias and R is the junction resistance (assuming 
zero impedance of the external circuit). Only a small 
fraction a of this power, i.e., aV 2 / R, is dissipated into 
the ionic degrees of freedom in the junction due to the 
electron-phonon coupling. The value of a need s to b e 
determined from a microscopic theory ( Todorovl Il998h . 
Since the spectrum of modes of the junction is typi- 
cally discrete, one expects a minimal bias (we call V c ) 
necessary to excite the lowest-energy phonon mode of 
the structure, and hence a(V < V c ) — 6 . Therefore 
we write P = Q(V — V c )aV 2 / R, where O is the step- 
function. Now, if the power P were not dissipated away 
from the junction, the latter would heat up substantially 
and eventually break down. Therefore, there must be a 
heat current Iq which dissipates this power into the elec- 
trodes. Since the leads are much bigger than the junction 
and are three-dimensional in nature, one can assume that 
this energy is carried away a t a bulk rate Iq = athTeS 



( Ashcroft and Merminl . Il976l ) , with T e s an average effec- 
tive temperature of the junction ions and o~th the lat- 
tice thermal conductance. At steady state the condition 
P = Iq then yields for the effective temperature 



Tcff = Q{V - V c ) 



CTthR 



1/4 



(13) 



Here, we have considered the bulk electrode temperature 
To = 0. If both electrodes are at finite temperature, then 
there is also a heat current ~ Ut^Tg flowing into the junc- 
tion, and hence the balance equation P — 2g ut — Iq gives 



6 In molecules, this bias may be very close to zero, due to the 
longitudinal "acoustic" mode of the whole molecule vibrating 
against the bulk electrodes. 
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V 



T cS = (T 4 +T4 )V* where T y = 0(V-V C ) (^) V4 
is the contribution to the temperature from the finite 
voltage bias. 

In the above considerations we have assumed that heat 
can be dissipated away from the junction rather eas- 
ily. The results may change depending on the heat- 
transport properties of the leads and the coupling be- 
tween the leads and the junction. For instance, if 
the leads are strongly disordered heat is carried away 
with a different expon e nt of the temperature difference 
( Yudson and KravtsovL l2003h . If the nano junction has 
poor thermal coupling to t he leads, or i n the presence 
of localized phonon modes ( Lepri etalX 120031 ). namely 
modes that have a very weak coupling with the bulk 
modes, then the local ionic temperature can reach very 
large values, even at relatively small biases. The reason 
is simple: in the above cases, due to the bias V, the 
current-carrying electrons are away from their ground 
state, and they are thus "seen" by the local modes 
of the nanostructure at an effective finite temperature. 
Thus, this situation provides the possibility for inelas- 
tic electron-ion scattering in an energy window ~ eV, 
with conse quent ion temperatures of the same order of 
magn itude (|Di Ventral . 120081: iTodorovl 11998k lYang et all . 
I2005D . We note that similar results were recently ob- 
taine d from microscopic considerations ( Mozvrskv et all . 
I2006T) . That being the voltage bias of 0.1 V would 

generate an effective temperature of ~ lOOOif. This 
seems to have been observed in atomic quantum point 
contacts at the breaking point ( Ward et aZ.L l2008h . Thus, 



good thermal coupling to the electrodes is essential for 
maintaining junction stability. 



3. Electron heating 



which was accompanied by a microscopic theory based 
on the quantum hydro dynamic equations for th e inter - 
acting electron liquid ( D'Agosta and Di Ventral . l2006h . 
is as follows. Assuming no electron-phonon interaction 
is present, to first approximation, the thermal electronic 
conductance of the electron liquid can be taken to be pro- 
portional to the temperature, a±h = jT e . Therefore, the 
heat current, given by Iq = cr t hT is quadratic in temper- 
ature, Iq — 7T e 2 . As in the case of local ionic heating, 
at steady state this thermal current has to balance the 
power dissipated in the junction, which is a small frac- 
tion of the total power of the circuit, P = aV 2 /R. One 
thus obtains 



T e = "f e - e V 



(14) 



where 7 e _ e is to be determined from a microscopic calcu- 
lation. Assuming the coefficient 7 e _ e weakly dependent 
on bias, this simple argument shows that the local elec- 
tron temperature grows linearly with bias. This result 
clearly hinges on the assumption that electronic heat is 
dissipated away from the junction at a bulk rate, which 
may not hold for all systems and under all experimental 
conditions. 

F rom a microscopic hyd rodynamic theory D'Agosta et 
al. ([D'Agosta et all . 120061 ) have also calculated the local 
temperature profile, T e (x), along the junction. From the 
maximal value of T e , an estimate of 7 e _ e was supplied for 
various junctions. For instance, for a 3D gold quantum 

o 

point contact (QPC) with effective cross section of 7 A", 
these authors evaluated j e - e (QPC) — 65 K/V. For a 
2DEG, assuming a cross section of 20 nm they found 
7 e _ e (2£>£G) = 1.2 x 10 2 K/V, suggesting that heating 
from inelastic electron-electron interactions is generally 
smaller than the corresponding heating due to electron- 
ion interaction. 



Up to now we have discussed the heating of the 
phonons in the junction due to their interaction with 
the current-carrying electrons. But what about the tem- 
perature of the electrons themselves? To be precise, 
we refer here to the temperature T e of the Fermi sea 
of electrons of the nano junction and those in its prox- 
imity. This temperature is affected by both inelastic 
electr on-electron interaction s and electron-phonon cou- 
pling ( D'Agosta et aU 12006). Clearly, the local electron 
temperature influences the local ionic temperature of the 
junction. However, accounting for both electron-electron 
and electron-phonon interactions is a challenging task. 
While attempts have been made to account for bo th in 
calculating charge transpor t dGalperin et al. . l2007bl) and 



recently even heat currents (|Liu et al. . 2009a), we are un 



aware of any calculation of the local temperature where 

these interactions are considered on equ al footing. 

In a recent work, D'Agosta et al (jD'Agosta et all . 
120061 ) have predicted the bias dependence of the local 
electron temperature in quasi-ballistic nanoscale junc- 
tions and its effect on ionic heating, treating the elec- 
tron liquid as a viscous fluid. The general argument, 



4. Ionic cooling 

A direct measurement of local electron temperatures, 
however, seems a very difficult task, and in fact we are 
not aware of such a direct method. On the other hand, lo- 
cal ionic temperatures are relatively easier to obtain (see 
Sec. IIII.Cp . It is then relevant to ask what is the effect 
of the local electronic temperature on the ionic heating. 
Since part of the total power dissipated in the junction 
goes into heating electrons via electron-electron interac- 
tions, that power is no longer available to induce ionic 
heating. Since the initial energy is always that of the 
current-carrying electrons, the ionic temperature must 
be smaller if electron heating takes place. The power of 
this electron-phonon scattering process can be assumed 
to have a form P e - P h — S( T" ft — T") with E a system- 
specific constant, and n > (jSchmidt et all . 12004V This 
ionic energy is then dissipated away from the junction. If 
we assume again a bulk dissipation law, Iq = a t hT^ s , for 
electronic temperatures much smaller than the ionic ones, 
the steady state condition P e - P h — Iq is satisfied by 
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E ~ o th and n = 4. By taking into acco unt a background 
temp erature T we then get the relation ([D'Agosta et all 
2006) 



T cff =(T 4 + 7p 4 



F 2 



7e 



,1/4 



(15) 



which is valid for V < {^e-ph/le-ej 1 '■ The meaning of 
Eq. [15] is that at sufficiently large biases, the effective 
phonon temperature is reduced, i.e., the phonons "cool 
down" . As we will discuss in the following Sec. IIII.Cl this 
result has been recently confirmed experimentally (see 
Figs. [7] and [8]). It is important to point out, however, 
that the exact power-law dependence in Eq. (|15[) and the 
value of the various coefficients may depend strongly on 
the details of the nanostructure and its contact with the 
leads. 

Another interesting idea to obtain reduced ionic heat- 
ing is to use a nanostructure with an appreciable Peltier 
coefficient. In this situation passing current through the 
junction would result in the cooling of one side of the 
junction, which may induce cooling of the molecule. The 
idea of local cooling of a junction has received renewed 
atten tion in recent years in the context of mes oscopic sys- 
tems (jGiazotto et aZ.l.l20 06; Sair a et all 20071) andmolec- 
ular junctions (iGalperin et all 2009bt iMcEnirv et all 
l2002HPistole^l2009tlZippilli et all 12003) . While the de- 
tails vary, the main concept is the same: the system is 
tuned in such a way that hot electrons (i.e., those with 
large kinetic energy) find it easier to tunnel through the 
junction, thus depleting the lead up-stream in voltage 
from hot electrons, thus cooling it. The cooling of the 
molecule is achieved either by its proximity to a cooler 
lead, or in more subtle cases, by the fact that elec- 
trons "borrow" energy from the localized phonon modes 
to assist transport, th us cooling them in the process 



(jGalperin et al 



p ort, tn 
., 2009b) 



C. Heating in current-carrying nanostructures: experiment 

Despite the difficulty in measuring directly local tem- 
peratures of nanoscale systems, we have witnessed much 
progress in this direction over the last years. The first 
concepts of local temperature meas urements are reviewed 
by Cahill et al. ( Cahill et aZ.I . l2003h . Especially notewor- 
thy are experiments where a thermocouple (serving as 
a thermostat) is mounted on top of an STM tip, thus 
creating a "scanning thermal microscope" (SThM) . This 
device was then used to study the local temperature of a 
carbon nanotube placed on a substrate (see Fig. EJ) - The 
authors of this work discuss several possible shortcomings 
and limitations of SThM studies: the dependence of the 
measured temperature on topography of the sample and 
surface chemistry, the fact that the tip itself might per- 
turb the sample (e.g., via near-field radiation, or by effec- 
tively cooling it), only surfaces can be measured, some of 
the heat is delivered through the air between the sample 
and tip, etc. These issues render this method hard for 
quantitative analysis, although some progress has been 
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FIG. 6 (Color online) Scanning thermal microscopy (SThM) 
of a lOnm diameter multiwall carbon nanotube. (a) Full ther- 
mal image, (b) A cut alon g the nanotube. (c ) A cut across 
the nanotube. Taken from i|Cahill et cZI. I2003I ) . 



achieved (|Grover et all \200& iKim et all l2008l ). We are 
unaware of any theoretical work (other than the one pre- 
sented in this review) which has been directly related to 
SThM measurements. 

In mesoscopic systems (e.g., quantum dots etched in 
2D electron systems) which are of typical sizes of mi- 
crons, tremendous advance in local thermometry has 
been achieved, as summarized in the thorough review by 
Giazotto et al. ( Giazotto et all [20061. In these systems, 



thermometry is achieved by analyzing some temperature- 
dependent function (current, conductance, etc.) from 
which, by using known properties of the electronic sur- 
roundings, the temperature can be extracted. To give 
a specific recent example, by analyzing the derivatives 
of the current a s a function of temper ature and voltage, 
Hoffmann et al. ( Hoffmann et aLl . l2009f) were able to mea- 
sure the temperature gradient across a current-carrying 
quantum dot of 15nm length, with the conclusion that 
the heat flow is mediated by phonons in the quantum 
dot. 

Other options for measuring the local temperature are 
available. One method is to study the force at which a 
molecular junction breaks as a function of current. The 
idea behind this method is that the higher the temper- 
ature of the structure, the less external force is need ed 
to break it (|Huang et all l2007t iTsutsui et all . l2008al lbh. 
From this force one can then extrac t an effective temper- 
ature. For example, Schulze et al. (Sc hulze et al. I. I2008D 
have studied the breakdown of a molecular junction com- 
posed of a Ceo molecule, and showed directly that better 
cooling of the junction is achieved when the coupling be- 
tween the molecule and the leads is improved. 

Using the above method, in a recent serie s of experi- 
ments, Huang et al. ([Huang et aLl . l2007l [20061 have stud- 
ied the local temperature of single-molecule (alkanethi- 
ols) junctions as a function of voltage bias. Their results, 
shown in Fig. [JJ (points) indicate that with increasing 
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FIG. 7 (Color online) Effective temperature of a molecular 
junction, for three different types of molecules n-alkanedithiol, 
with n = 6 (squares), n = 8 (circles) and n = 10 (triangles). 
The solid lin es are theoret i cal es timates from Eq. 1151 From 
Huang et al. (|Huang et all 120071 ). 



voltage, the local temperature first increases, saturates, 
and then slightly decreases. This is in agreement with 
the prediction of Eq. [15] (solid lines) , and suggests that 
electron-electron interactions indeed occur in these junc- 
tions. The same experiment also confirms that longer 
alkanethiol molecules heat up less due to increased el ec- 
tronic resistance, at fixed voltage ( Chen et~al\ . l20 05b). 



ature has been suggested recently. It makes use of 
Raman spectroscopy and it w as first applied to the 
study of a suspended n anotube ( Bushmaker et q?J . l2007t 
iDeshpande et all [2009) . In these measurements, the lo- 
cal temperature was deduced from the shifts in the local 
Raman G+ and G_ bands of the nanotubes. The au- 
thors compared two nanotubes of lengths 2/im and 5/im, 
and found that the longer nanotube was less heated, an 
effect which was attributed to the thermalization of hot 
phonons at the center of the nanotube. 

Along the same lines, the local temperature of a molec- 
ular junc tion has b e en in vestigated via Raman spec- 
troscopy ( Ioffe et all l2008f ). The idea is that the ra- 



tio between the Stokes and anti-Stokes intensities is di- 
rectly related to their non-equilibrium populations in the 
presence of electronic current. The method h a s been 
discussed theore tically in detail ( Galperin et all l2009at 
Ioffe et all 120081) . In Fig. [5] the effective temperature is 
plotted as a function of voltage bias, for different Raman 
modes. From the figure one can see that although the 
temperature is slightly different between different modes 
(due to the different electron-phonon coupling strengths 
and symmetries), the overall voltage dependence shows 
roughly similar features. Note also the apparent decrease 
in the local temp erature, which is in line with the results 
of Huang et al. ( Huang et aZ.I . l2007h . This study, along 
with the one described before, indicate that local Raman 
spectroscopy may serve as a valuable tool for the study 
of local temperatures at the nanoscale. 
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FIG. 8 (Color online) Effective temperature of a molecu- 
lar j unction measured using Raman scattering; from Ioffe et 
al. (jloffe et all l2008h . The different points correspond to 
different modes, and while the temperature is slightly dif- 
ferent, the overall voltage dependence shows roughly similar 
features. Note also the apparent decrease in the local tem- 
perat ure with bias, which is in line with the results of Huang 
et al. (|Huang et all 120071 ). 

An alternative method to study the local temper- 



D. Fourier's law at the nanoscale 

Let us now end this section with a somewhat differ- 
ent issue, that of the onset of Fourier's law, Eq. (|T|), in 
nanostructures. As previously noted, in the context of 
nanoscale junctions, there is not much point in discussing 
the scaling of the thermal conductivity n, which pertains 
to an asymptotic relation, strictly valid in the limit of 
large system lengths. Therefore, here we limit our dis- 
cussion to the simple question of what is the temperature 
profile along the junction. 

In the context of Fourier's law, this question has been 
widely studied for both classical (iLepri et al. . 2003) and 
quantum systems (jMichel et al ., 2006). The main focus 
has been on e ither spin-chains (i-e . . Isin g-like models) 
(jMichel et all 120061: IWu and Segal . 120081) or harmonic 
oscillator chains. The local temperatures are usually 
evaluated by c alculating the averages of some local en- 
ergy operato r s dMeiia-Monasterio and Wichterichl . l2007t 
iMichel et all 120061 120031: ISaitol. l2003l). o r by using self- 
consistent thermal baths dDharl . l2008l iDhar and Rovl . 
120061 : iJacauetl 120091: iRovl . 120081) . In the first case, one 
assumes that the local energy is re lated to the tempera- 
ture vi a a local Boltzmann relation ( Dubi and Di Ventral . 
l2009bl ), or directly propo rtional to the temp erature via 
a local equipartition law ( Michel et all [2003) . The dis- 
advantage of this method is two-fold: (i) there is some 
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arbitrariness in choosing the local energy operator, since 
one can represent the same Hamiltonian in different ways 
( Wu and Segall . l2009h . and (ii) this method assumes from 
the outset that the system is in a local thermodynamic 
equilibrium, which may not always be the case. 

In the second approach, the system is attached to local 
heat baths. The heat current between the junction (or 
quantum wire) and the local baths is calculated, and the 
temperatures of the heat baths are determined in such 
a way that the heat current between the wire and the 
baths vanishes. This method was recently described in 
detail and applied to a quantum chain of non-interacting 
harmonic oscillators (iRovl . 120081) and a ch ain o f quan - 
tum dots (|Jacauetl l2009h . For instance, in (jRovl . 120081) . 
using quantum Langevin equations, the local tempera- 
ture as a function of position was calculated for differ- 
ent chain lengths and for different coupling between the 
wire and the local baths. The conclusion of this work 
is that the coupling between the wire and the baths de- 
termines a length scale (mean- free path), and the heat 
transport crosses over from a diffusive regime (uniform 
temperature gradient) to a ballistic regime (uniform tem- 
perature, vanishing gradient) depending on the system 
length being longer or shorter than the mean free path, 
respectively. Since the dynamics of the system is calcu- 
lated in the presence of the local baths, this shows that 
the properties calculated (i.e., local temperature) pertain 
to the combined system of quantum chain and thermal 
baths, and thus naturally depends on, e.g., the coupling 
strength between them. 

Recently, a method has been suggested to calculate 
the l ocal temperature of electron s in a nanoscale junc- 
tion ( Dubi and Di Ventral 12009(31(1 ) . Its starting point is 
the stochastic Schrodinger equation (see Eq. ([2"B]0 . which 
for non-inter acting electron s redu ces to a quantum mas- 
ter equation ( Pershin et all 120081) . In this approach the 
finite electronic system is coupled to two local heat baths 
at the edges of the system, in similarity to the study pre- 
sented above for a chain of harmonic oscillators. In order 
to evaluate the local temperature, the definition we in- 
troduced in Sec. III. Al has been used. Namely, a third 
environment is coupled locally to the system at the posi- 
tion where the temperature needs to be evaluated. The 
properties of the system are then evaluated twice: once 
with the additional environment (so-called "tip", as it 
mimics, e.g., the operation of a thermostat mounted on 
an STM tip) and once without the probe. The tem- 
perature of the probe is then varied (floated) such that a 
minimal change in some local (or global) properties of the 
system, such as its local electron density, occurs. A scan 
of the local temperature of the whole system can then 
be obtained with this method. The advantage of this 
approach is that it can, in principle, be implemented ex- 
perimentally, and it provides the local temperature of the 
electrons without further scattering from other sources. 
In addition, it can be shown analytically that the above 
definition reduces to the standard thermodynamic tem- 
perature in limiting cases, for instance, in local equilib- 



rium (see also ( Di Ventra and Dubil . l2009f) ) or for two- 
level systems. 

For the case of a wire coupled to two electrodes at dif- 
ferent temperatures, it was found that the local temper- 
ature of the wire may exhibit qua ntum oscillations for 
interm ediate lead- wire couplings ( Dubi and Di Ventral 
l2009dl) . Similar oscillations were later observed for a 
driven quantum wire (jCaso et and reflect the 

quantum coherent nature of the system. When the lead- 
wire coupling is large enough a uniform temperature en- 
sues. In this limiting (ballistic) case, one also finds that 
the non-equilibrium distribution function of the system 
is an average of the distribution functions of the left and 
right baths. The fact that the temperature is uniform in 
the wire demonstrates the known result that for a clean 
system, Fourier's law is invalid. 

In order to reconstruct Fourier's law (with an associ- 
ated temperature gradient), diagonal disorder was intro- 
duced in the wire (which localizes the electronic wave- 
functions), and the lo cal temperature was averaged over 
disorder realizations (Dubi and Di Ventr 1 l2009dY It 
was found that for large enough disorder, a local uni- 
form temperature gradient ensues, giving rise to Fourier's 
law. This result was interpreted in terms of an ef- 
fective thermal length which controls the scale o f the 
temperature gradient ( Dubi and Di Ventral . l2009bl ) . By 
adding the effect of dephasing th e model wa s also able 
to explain the results by Roy (|Rovl H008) described 
above. We finally conclude that for the above model 
the on set of Fourier's law coincid es with the onset of 
chaos ( Dubi and Di Ventral l2009cD . This has also been 
found in oth e r mod el systems ( Gaul and Buttnerl 120071 : 
Michel et all |2006f). but not in all cases (jLepri et all 
Thus, this result does not appear 



1997; Li et al\ 



2002). 



to be universal. 



IV. THERMOPOWER 



A. Introduction and basic definitions 



In this section we discuss the concept of thermopower 
in nanoscale junctions. As prototypical examples that 
show all main features of the problem we will be fo- 
cusing mainly on experiments in molecular junctions 
( Reddv et all l2007f ) and briefly mention experiments in 
mesoscopic systems and nanowires. The thermopower 
phenomenon corresponds to the case in which a temper- 
ature difference at two sides of a given junction induces 
a voltage drop across it. From a technological point of 
view, this effect is of great importance, since it may be 
used to recover part of the heat wasted in physical pro- 
cesses and generate electrical power with no moving me- 
chanical parts. It is also of basic scientific interest, since, 
by combining energy and charge flow, it may encode in- 
formation about the system dynamic s which is un avail- 
able in charge transport experiments (Segal [2005). 

The configuration we have in mind is again a junction 
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composed of two leads separated by a nanoscale element 
- a quantum dot, a molecule, nanotube, etc. Consider 
such a junction, where the two leads are held at different 
temperatures, Tl and Tr. The corresponding tempera- 
ture difference AT = Tr — Tl gives rise to both a heat 
current (discussed in Sec. UTJ and a charge current. If the 
circuit is closed, after a transient time charges accumu- 
late on one side of the junction and deplete on the other, 
so that a zero charge current is achieved and a voltage 
drop across the junction is formed. If the circuit is open 
(namely it is connected to an electron source), and a volt- 
age difference AV is applied between the two leads with 
appropriate sign, a bias-induced electric current can can- 
cel out the thermally-induced current. Note, however, 
that given a temperature difference the two procedures 
may not yield the same voltage difference. In fact, the 
voltage difference may also depend on the location along 
the system where it is probed. 

The thermopower S is defined as (minus) the amount 
of voltage AV at the state of vanishing current, 



S = 



AV 
AT 



(16) 



1=0 



in the limit of AT -> 0. 

This definition can also be understood from the current 
expressed in linear response. This is defined as 



I = GAV + L T AT, 



(17) 



where G is the electrical conductance and Lt is a re- 
sponse coefficient related to the energy flow. From this 
expression one readily sees that S = Lt/G. Therefore, 
in order to determine S, one has to determine the con- 
ductance G and the thermal response Lt- 

Before we proceed to discuss different theoretical and 
experimental aspects of the thermopower, it is impor- 
tant to ask the following question: is knowledge of the 
thermopower S sufficient to design devices that operate 
as efficient heat- voltage converters? In fact, in a real 
device, it is not at all clear that the system is, under 
the given experimental conditions, in the linear response 
regime. Nor it is obvious th at the best conversion should 
be achieved in that r egime (|Dubi and Pi Ventral 12009^ ; 
lEsposito et"aZl . l2009h . 

To this end, it is useful to define the unit-less "figure 
of merit" , 



ZT = 



GS 2 
<Jth/T 



(18) 



where T is the temperature of the system 
(|Mahan and Sofol . fl996l ) . The quantity ZT describes the 
efficiency of a real device or material as a thermo-electric 
converter. While an exact relation between ZT and 
thermodynamic efficiency is available ( Mullen , |2008| ) . this 
choice can be intuitively understood: S measures how 
large a voltage drop can develop for a given temperature 
gradient, G measures how easy charges can cross the 



junction to generate that voltage drop, and <ith measures 
how hard it is to maintain a temperature gradient. 

It is commonly stated that for applications, one must 
achieve ZT » 1 (in fact ZT > 4 would already be 
a great advance (iMaiumdari . I2004D ). However, such a 
situation is hard to obtain: in most cases the electrical 
conductivity a and thermal conductivity k are related vi a 
the Wiedemann- Franz law (jAshcroft and Mermmlll976T) . 
which states that 



K 

a 



3e 2 



T, 



(19) 



with the quantity in parenthesis commonly referred 
to as the "Lorenz number" . This means that it 
is difficult to increase a and S without also in- 
creasing k, and vice versa. However, deviations 
from the Wiedemann-F ranz law have been observed 
( Applevard et all l2000h and discussed the oretically in 
vario us systems, includin g Luttin g er liquids dGarg et al 



20091 iKane and Fisherl Il996t iKubala et all I200S : 



Li and Orignad l2002t iMurphv et all 120081 iReiec et 
20021) . These deviations are attributed to interac- 
tions, where the simple single-particle theory fails (see 
Sec. HV.Cl on theoretical methods), and are exactly what 
is required in order to increase the efficiency of thermo- 
electric devices. 



B. Experiments on thermopower at the nanoscale 

Measurements of thermopower are conceptually easier 
than those of thermal conductance: one applies a tem- 
perature gradient across the junction, and measures the 
ensuing voltage in a closed circuit when the transient 
current vanishes. Or, in an open circuit, one supplies a 
voltage to compensate for the thermally-induced current. 
The slope of the resulting voltage-temperature gradient 
curve gives the thermopower. However, in an actual ex- 
periment, particular care needs to be taken to extract this 
quantity. The reason is because the voltage probe that 
is connected to the system in order to measure the ther- 
mopower is necessarily invasive, since the applied thermal 
gradient would induce, locally at the voltage probe con- 
tact, an extra voltage difference. This extra effect needs 
to be subtracted to get the actual thermopower of the 
nanojunction. In addition, the ensuing voltage (includ- 
ing its sign) is very sensitive to the junction geometry 
and thus may fluctuate considerably in an actual exper- 
iment, providing non-trivial distributions of the voltage 
as a function of thermal gradient, from which a single 
voltage value may not always be easy to extract. 

Before reviewing some recent experiments on nanoscale 
junctions, it is of interest to briefly survey some of the 
older experiments on mesoscopic systems as well. We 
point out that while most of these experimental results 
may be understood in terms of a linear response scatter- 
ing theory (see next Sec. HV.Cj) . some recent results, such 
as the appearance of additional peaks in the distribution 
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of induced voltages versus temperature gradient (Scheib- 
ner et al, 2005), are yet to be completely accounted 
for. The discovery of pronounced mesoscopic effects 
such as Coulomb blockade and conductance quantization 
prompted the study of thermopowe r in quantum point 
conta c ts i|van Houten et all Il992bt Molenkamp et al. 



Godii n et al. 



1999: 



1994 [l990D and quantum dots 

Scheibner et all 120051 : IStaring et all Il993l ). These de- 
vices are defined by depositing gates on top of a two- 
dimensional electron gas formed in a semiconductor in- 
terface (typically GaAs/AlGaAs). Heating of one side of 
the device is achieved by passing current through it with 
consequent Joule heating and temperature rise. Most 
of the results of these experiments are well understood 
within the simple, s ingle-p article picture of thermopower 
( van Houten et al .1. Il992bl) . which will be described be- 
low. 

Another, more recent batch of thermopower experi- 
ments are those conducted on nanowires, namely wires 
with nanoscale diameter, but extending in the lon- 
gitudinal direction as long as a few microns. Var- 
ious experime nts were performed on wires of differ- 
ent materials dBoukai et al. . 20061 20071 Duarte et al. 



2009 1 : iHochbaum et aZ.I . l2007i iLee eVaLl.l2009HSeol et al. 



2007), as well as carbon nanotubes (IKong et all 12005 
Small et all . [2003: Sumanasckcr a et «Ll . l2002l ). These ex 



periments suggest that in these systems it is possible to 
either increase S (by designing the system to have an 
increased electronic density of states) or reduce the ther- 
mal conductance independently by, e.g., designing a sys- 
tem which is smaller than the phonon mcan-free-path 
but still larger than the corresponding mean-free-path 
of the electrons or holes, thus increasing the figure of 
merit. Specifically, in Si nanowires these are obtained by 
the combined effect of the change in phonon spectra and 
enhanced scattering off the boundary, both having little 
effect on the e lectronic part (see recent review by Rurali, 
( Ruralil . l2010h ). Along similar lines, boundary effects 
seem to highly reduce the thermal conductan ce but leave 
the c harge conductance roughly unchanged ( Maiumdarl . 
I2004D . 

A set of experiments which is of great interest both 
from an academic and technological point of view, are 
those performed on j unctions of nanometer length, s uch 
as metallic contacts ( Ludoy ji and Ruitenbeekl. 1 19 99) or 
molecul ar junctions llBaheti et all 2008: Malen et all 
l2009aH bl: lReddv et all [20071: iTan et aZ.l . l2010l) . The latter 
ones are of interest since, as we will discuss in Sec. IIV.C1 
theoretical arguments suggest that molecular junctions 
may exhibit large thermopower. In these latter experi- 
ments, a gold STM tip is placed on top of a gold sub- 
strate which is covered with various molecules. As the 
STM tip touches (and is attached to) a molecule, a ther- 
mal gradient is applied and the thermopower is measured 
by applying a voltage so that no current passes through 
the junction (see upper panel of Fig. [S]). This proce- 
dure is repeated many times and a histogram of the volt- 
age required to achieve vanishing current is obtained (for 
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FIG. 9 (Color online) Upper panel: schematic representation 
of thermopower experiments on molecular junctions, (a-c) 
Distribution of thermo-voltages obtained at different temper- 
ature gradients. Note the widening of the distributions and 
their non-trivial structure, (d) The most-probable thermo- 
voltage obtained from (a-c) as a function of the temperature 
gradient. The derivative of the linear fit of this curve yields 
the thermopower S. (e) S of various molecules, i n terms of the 
mole cule's length (adapted from Reddy et ai. dReddv et all 
[2003)). 



different temperature gradients, AT = 10K, 20K, 30K, 
Fig. Ufa-c)). The authors of this work then plot the volt- 
age with maximum probability (i.e., the position of the 
peak in the voltage histograms) as a function of AT, and 
by fitting the resulting curve with a linear fit the ther- 
mopower is obtained (Fig.[^d)). These experiments were 
performed with various kinds of molecules, and interest- 
ing phenomena such as a length dependence of the ther- 
mopower (Fig. in^e)) or strong dependence of the molecu- 
lar end-groups were observed. The experiments indicate 
that indeed molecular junctions have favorable thermo- 
electric properties, suggesting that devices incorporating 
molecular junctions or arrays (for instance metallic plates 
separated by a molecular layer) may be good candidates 
for thermoelectric applications. 

Let us, however, point out some features of the exper- 
iments which at present do not have a satisfactory theo- 
retical explanation. For instance, as seen in Fig. Htc), 
the voltage histograms have a well-defined structure, 
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with a not so negligible secondary peak near AV = 0. 
Note also that the distributions cross the AV=0 line 
into negative values, not shown in the figure. An addi- 
tional feature of the histograms is their apparent widen- 
ing with increasing temperature gradient. These fluc- 
tuatio ns effects have been recently studied experimen- 
tally (|Malen et all l2009ah and are attributed to varia- 
tions in contact geometry and orbital hybridization, as 
well as intermolecular interactions, in accord with theo- 



retical studies (jDubi and Pi Ventral [2009dT) . 

The analysis of the above results has been done based 
on the single-particle (non-interacting) Landauer ap- 
proach to thcrmopower. As it will be discussed in 
Sec. IIV.C1 in the linear response single-particle the- 
ory of thermopower, S can be simply related to the 
electronic properties of the junction, and specifically 
in molecular junctions, to the position of the elec- 
trochemical potentials of the leads with respect to 
the gap between the highest occupied molecular or- 
bital (HOMO) and the lowest unoccupied molecular or- 
bital (LUMO). Since the position of the HOMO-LUMO 
gap affects charge trans port in molecular junctions 
(jNitzan and Ratnerl l2003h . measurements of S were sug- 
gest ed as a way to probe the energy position of thes e lev- 
els (|Baheti et all 120081: iPaulsson and Dattal. 120031). I n 
the experiments (|Baheti et aLl . l2008HReddv et aLl . l2007l) . 
a comparison of the thermopower and conductance with 
numerical simulations using ground-state DFT within 
the Landauer approach has been performed. From this 
comparison it was then concluded that the position of the 
HOMO-LUMO gap can be thus determined. This anal- 
ysis, however, raises several questions. The applicability 
of a linear-response single- particle theory was questioned 
by the authors themselves ( Reddv et aZ.1 . 12007^ . following 
the fact that the temperature differences can be tens of 
degrees Kelvin. In fact, one could ask if this is the small- 
est energy scale in the experiment. Specifically, is this 
energy smaller than, say, the coupling energy between 
the molecule and the substrate? The answer to this ques- 
tion is unclear, especially in light of the large error bars 
shown in Fig.[9jd). In addition, the non-trivial structure 
of the fluctuations in the voltage histogram implies that 
non-equilibrium effects may come into play, which are 
not taken into account in the linear response theory. Fi- 
nally, electron-electron and electron-phonon interactions 
may play a crucial role in this problem. 

Despite these open questions, the experiments de- 
scribed above are very impressive and important for the 
field and there are many interesting future directions to 
which they could be taken. For example, it would be 
interesting to study the change in the width of the distri- 
butions and their structures as the overall temperature 
is reduced. This would determine, e.g., if these distri- 
butions are due to static or dynamic effects. Another 
interesting direction would be to study, for several molec- 
ular structures, not just the most-probable voltage, but 
the real (statistical) average of the distributions, and in- 
fer from this whether the resulting thermopower displays 



the same features as reported above (e.g., length depen- 
dence, etc.), and whether this quantity matches calcula- 
tions based on single-particle theories. 



C. Theoretical methods 

In this subsection we describe the present theoreti- 
cal methods available to describe the phenomenon of 
thermopower. The most common one is based on 
the Landauer approach with its most common imple- 
mentation within ground-state DFT. As we will dis- 
cuss, this approach has several advantages, being rather 
computationally straightforward, and having a rather 
simple physical interpretation. However, we will ar- 
gue that in many cases of actual experimental inter- 
est, it may be inadequate, since it is based on the 
usual assumptions of scattering theory of non-interacting 
electrons. In addition, as we have emphasized also 
in Sec. III. Oil the use of ground-state DFT is ques- 
tionable in an intrinsically non-equilibrium problem as 
that discussed here. We will then introduce an ap- 
proach based on the theory of open quantum systems, 
which is ideally suited for the present problem and 
can, in principle, account for interactions (beyond mean- 
field). The latter point has its most practical imple- 
mentation in an exten sion of time-dependent DFT to 
open quantum systems (Ip'Agosta and Pi Ventral . l2008bl 
iDi Ventra and D'Agostal . l2007fl 



1. Single-particle theory of thermopower 

The starting point for calculating the thermopower 
within a single-particle pictur e is the Landa uer expres- 
sion for the electrical current (jButchen . fl99(l . 



/ = 



deT(e) [/ L (e) - / R (e)] 



(20) 



where T(e) is the transmission coefficient at energy e and 
fh R are the Fermi distributions of the left and right leads. 
In the limit of small bias and temperature gradient (i.e., 
| AT | << T and \eAV\ « ^, where T is the background 
temperature and \x is the equilibrium chemical poten- 
tial ) the distributi on functions are approximately given 
by (lButched . fl990h [i = L,R) 



f(e,fj,,T)± — (n-Hi)T-^{e-n) 



(21) 

where now /(e) is the equilibrium distribution, and the + 
and — sign correspond to which electrochemical potential 
is higher or lower in energy with respect to the equilib- 
rium chemical potential. Inserting this into Eq. (|20[) and 
equating the current to zero, one obtains 

1 r oo d £ T(e)( e -^)(-,r( e )) 
MT) "eT r oo deT(e)(-/'(e)) " ^ 
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Already from this result several features may be seen. 
First, since at T = we have —/'(e) = 8(e — /i), the 
numerator of S vanishes and S(T = 0) = 0. Second, 
even at finite temperatures /'(e) is symmetric around /i, 
and therefore S = unless 7~(e) is not symmetric around 
(i. This is similar to the condition in bulk materials that 
requires the particle -hole symmetry be broken to have a 
finite thermopower (jAshcroft and Merminl . 119761) . 

One can further simplify S(T) by taking the low- 
temperature limit and by assuming that there are no 
electronic resonances close to the equilibrium chemical 
potential. Using the So mmerfeld expansion to first o rder 
around fi(T = 0) = e_p ( Ashcroft and Merminl . [l976h one 
has 



Wiedmann- Fr anz law, thus r educin g the agreeme n t with 
experiments dKubala et all . 2008 ; Limde et all 2006 ; 



T(e)(e- M )/'(e) 



d 2 [T(e)(e-/i)] 



de 2 



y4^T'(e) 



(23) 



(where the second derivative comes from integration by 
parts) and one arrives at the expression for the ther- 
mopower, 



^k B dj<r(e)) 
S = k B T 

6 e de 



(24) 



which is similar to Mott's semiclassical fo rmula 
for bulk metals (lAshcroft and Mermi 3, fl97l 
iLunde and Flensber 1 l2005h . We stress once more 
that this approximation is only valid at low tempera- 
tures and away from transmission resonances, so that 
the variation in T(e) is small. 

The advantages of using the Landauer formalism 
are evident: it provides both a simple interpretation 
of the thermopower in terms of single-particle prop- 
erties such as the transmission coefficient T(e), and 
a rather straightforward computational procedure. In 
fact, one only needs to determine the transmission 
coefficient T(e), which can be done as discussed in 
Sec. III.C.ll These reasons have made this approach 
extremely popular and widely used. An early use of 
Eq. (I24[) is in the study of the r mopow e r in quantum point 
contacts (Ivan Houten et al\ , Il992b ; Molenkamp et all 
I1994L Il990h and quantum dots (jStaring et all Il993h. 
In these mesoscopic systems, a gate voltage is used 
to tune either the width of the quantum point con- 
tacts or the energy levels in the quantum dots, giv- 
ing rise to quantized conductance and Coulomb block- 
ade. It turns out that in the cases above, the Lan- 
dauer approach yields reason ably good agreement be - 
tween theory and experiment ( Molenkamp et a.7] . ll990T> . 
and knowledge of T(e) reasonably describes both the 
conductance and the thermopower. This would naively 
suggest that for these types of systems, the above 
single-particle picture accounts for most of the ther- 
mopower. However, more recent investigations which 
include effects of interactions, show that in both types 
of systems interactions may induce deviations from the 



Turek and Matveevl 12001 iTurek et all 120051 : iZhangl 
20071 ). suggesting that the agreement in the single- 
particle case may be the result of cancelation of errors. 

In fact, despite its simplicity, the above approach 
suffers several shortcomings of particular relevance in 
nanoscale systems. The most prominent is the fact 
that it is formulated for non-interacting electrons. This 
means that any inclusion of interaction effects di- 
rectly iri^JTXlL-SSiL^HliL.. ^—^ , 116 a t the mean-field 
level ( Vignale and Di Ventral l2009t ) . To correct this, one 
should abandon the Landauer formula for the current, 
and, alternatively, use expressions for the currents ob- 
taine d by using, e.g., the NEGF (IMeir and Wingreenl . 
Il992h or rate equations (|Koch et all \2004 . To our 
knowledge, in its fully interacting form NEGF was never 
employed to study the effects of electron interactions on 
the thermopower of molecular junctions. 

Another limitation of the Landauer approach is the 
erroneous result it supplies in the zero-coupling limit. 
To demonstrate this, consider the simplest model for a 
nanoscale junction: a single resonant level symmetrically 
coupled to leads with spinless electrons (adding spin sim- 
ply introduces a factor of two) . The transmission is given 
by a Breit-Wigner formula, T(e) = F 2 /(r 2 + (e - e F ) 2 ), 
where T is the lead-induced level broadening. Substitut- 
ing into the expression for S (Eq. (|M|) ) and taking the 

limit of T — > gives a finite value, S = -2^-% — ^— . 
However, if one would consider a real device, it is 
clear that by detaching the leads would result in no 
temperature-induced voltage drop. The reason for this 
discrepancy is simple: in the linear response calculation 
one assumes that the temperature difference is the small- 
est energy scale, yet in the limit T — > 0, T becomes com- 
parable to such a scale, and the approximation breaks 
down. One should thus be careful both in using per- 
turbation theory in the coupling between the leads and, 
say, a molecule in the junction and in comparing such 
calculations to experiments (see Sec. IIV.C.2I) . 

Much of the recent theoretical work on thermopower 
has been devoted to molecular junctions. Before we 
review some recent results, it is important to under- 
stand the origin of the specific interest in such sys- 
tems, which may be understood from analyzing the Lan- 
dauer formula (|24[) . In a molecular junction, the Fermi 
energy of the leads is placed somewhere between the 
HOMO and the LUMP (i e., in the HOMO-LUMO gap) 
()Nitzan and Ratnerl . 120031) . The question is where ex- 
actly? The answer to this question cannot be answered 
by studying the conductance (or transmission) alone, 
which can be demonstrate d through a simple example 
(jPaulsson and Dattal . |2003). Consider such a molecular 
junction, with HOMO and LUMO energies enoMO an d 
clumOi respectively. The transmission function can be 
modeled by a double Lorenzian, corresponding to tunnel- 
ing via each of these levels independently, and assumes 
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the following form, 
T{e) = , rLFR 



r L r R 



r 2 + (e F - ehomo) 2 T 2 + (e_F - elumo) 2 

, (25) 

where Fl,r is the level broadening due to the left (right) 
lead and T — (IY + Tr)/2. (For simplicity, we as- 
sume it to be the same for the two levels.) The re- 

suiting thermopower (in units of ^--^-T), along with 
the transmission coefficient, is plotted in Fig. [10] (tak- 
ing r/(e LU MO - fHOMo) = 0.1). As seen, according to 
this simple model, for a given value of transmission be- 
tween elumo and £homo, there are two values of the 
Fermi energy which provide a solution to Eq. (|25|) . and 
hence conductance alone does not suffice to determine the 
position of the Fermi energy. From the same model, how- 
ever, one would infer that the sign of the thermopower 
is determined by the position of the Fermi energy from 
the center of the HOMO-LUMO gap, similarly to the 
fact that the sign of the thermopower in bulk materials 
is determined by whet her the conductance is domi nated 
by electrons or holes ([Ashcroft and Merminl . Il976j) , and 
therefore the sign of thermopower distinguishes between 
the two Fermi energies which solve Eq. (|2"5"|) . 
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FIG. 10 (Color online) Transmission (solid line) and normal- 
ized thermopower (dashed line) as a function of the position 
of the Fermi energy with respect to the HOMO-LUMO gap, 
based on the Landauer formula, Eq. (|24|) . with approximation 
Eq. (25]). 

This idea, along with the prospect of using molecular 
junctions as efficient thermoelectric devices, has gener- 
ated much theoret ical interest. To give a few examples, 
Segal (Segal 120051 ) showed that by measuring the ther- 
mopower one can distinguish between different electron 
transport mechanism s. Thermal and v ibrational effects 
were studied in detail ( Koch et Ml, [2004) using rate equa- 
tions, and it was shown that at low temperatures the sig- 
nature of the vibrational mod es on the thermopower can 
be measured. Murphy et al. ( Murphy and Moorel |2007| ) 
used rate equations to study the optimization of the fig- 
ure of merit of a molecular junction, in similarity to the 
optimization of the figur e of merit in bulk thermoelectrics 
(jMahan and Sofd . Il996h . 



Much recent attention has been devoted to study- 
ing thermopower of molecular junctions using ground- 
state DFT to calculate T(e) combined with the 
Landauer formula (1221) (or its even more simpli- 



fied v ersion Eq. (l24l)1 dFinch et all [2009; 



2009; iLiu and Chen. 200 



2008: iPaulsson and Datt 



Liu et al 
2003k 



Ke et 



2009b; Miiller 



Pauly et all 



2008; 



Wang et all 120051; IZheng et all 120041 " In some cases, it 



has also been reported an impressive ag reement betwee n 
the theoretical results and experiments ( Ke et adl2009h . 



However, caution has to be applied in making such 
claims. In fact, if the system is away from linear 
response - and many experiments so far likely corre- 
spond to such case - given a temperature difference, set- 
ting I = in Eq. [5D] does not necessarily provide a 
unique solution for the potential difference. In other 
words, there may be more than one potential differ- 
ence AV that gives rise to the same AT (and hence 
several values of thermopower for the same tempera- 
ture difference), when / = 0, as it is permitted by 
the non-linearity of Eq. [20] In addition, as already 
emphasized previously, even if the single-particle equa- 
tions (f2"2"|) and (f2~4"|) were good starting points to describe 
the problem at hand, ground-state DFT is fundamen- 
tally flawed in the present context (even if one knew the 
exact ground-state functional) due to t he fundamental 
non-equilibrium n ature of the problem (iBushong et al 



non-equilibrium n atur< 
2005 : iDi Ventral [200 



Vignale and Di Vcntra 



ng 

.Di Ventra and Todorovl 
20091) 



2004; 



In this respect, even the 
interpretation of ground-state Kohn-Sham orbitals is 
questionable, since the latter ones are auxiliary quan- 
tities whose only role is to provide the correct density of 
the corresponding many-body system in its ground-state. 

To these limitations we must also add few more physi- 
cal issues. When a thermal gradient is applied to a junc- 
tion, the transient dynamics is fundamental in establish- 
ing the voltage difference that enters the definition of 
thermopower. Since the dynamical formation of local re- 
sistivity dipoles creates strong local fields at the junction 
(especially at the nanoscale), these fields influence the 
electron motion in a non-trivial way, and thus influence 
the long-time behavior of the carrier dynamics, even in 
the dc limit. T his is particularl y important away from 
linear response (IDi Ventral fcOQSL which may be the ex- 
perimental case. 

It is the self-consistent formation of these fields 
that makes the thermopower very sensitive (both in 
magnitude and sign) to atomic details, and thus 
to the contact geometry between the nanostructure 
and bulk electrodes, as demonstrat ed also experimen- 



tally ( Ludoph and Ruitenbe ek. 1999). This precludes an 
easy interpretation in terms of "electron" or "hole" exci- 
tations as in bulk metals, and thus an easy relation with 
single-particle states (such as the HOMO and LUMO) as 
the Landauer equation (|22l) would imply. All this points 
to the fact that, since the system is in dynamical inter- 
action with two different baths, one needs to go beyond 
the approximations underlying Eq. (I22[) . and consider an 
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open quantum system approach. 



2. An open quantum system approach 

The present authors have precisely explored the 
problem of thermo power within the t h eory o f open 
quantum systems ( Dubi and Di Ventral . l2009dl ). In 
analogy with the idea that ele ctrical current s may 
be studied using finite syste ms ( Bushong et all 



iDi Ventra and TodorovLl2004[ ). one can study a finite sys- 
tem in contact with two heat baths held at different tem- 
peratures (i.e., finite leads connected by a nanoscale con- 
striction, either a molecule, wire, etc.). If the system 
has a finite thermoelectric response then charges would 
flow between the leads until the ensuing electric potential 
"balances" the thermal gradient, and a charge imbalance 
is created between the two leads (which is related to the 
thermo- voltage via the Poisson equation). Note that in 
this approach the system is allowed to find its own charge 
distribution via the transient dynamics (unlike a static 
approach where a static distribution is imposed a priori 
via boundary conditions) , and even when the charge cur- 
rent is zero an energy current is still present, as in the 
actual experiments. Then, by calculating the thermally- 
induced charge imbalance one obtains information on the 
thermoelectric response of the junction via the usual def- 
initions. This approach is also not limited to linear re- 
sponse thus providing information on the thermo- voltage 
even when the temperature gradient is not the smallest 
energy scale. 

An implementation of such an open system approach 
can also be formulated within time-dependent DFT, 
thus allowing to include, in principle, all possible dy- 
namical many-body effects in the thermopower (re- 
call in fact that given the baths that set the tem- 
perature differences, the ensuing electrostatic voltage 
is a well-defined functional of the density). Indeed, 
Di Ventra and D'Agosta have recently proved that 
if the bath-electron interactions are treated within a 
memory-l ess approxima t ion ( the thermal baths being 
Ohmic) (|Van Kampenl . 12001 ). then there is a one- 
to-one correspondence between the exact ensemble- 
averaged current density and the external vector poten- 
tial, therefore extending the theorem (and Kohn-Sham 
scheme) of time-depen dent current-DFT (TDCDFT) to 
open quantum systems (TD'Agosta and Di Ventral . l2008bt 
IDi Ventra and D'Agostal 120071 ). The framework for this 
theory (named Stochastic TDCDFT) is the stochastic 
Schrodinger equation, which describes a Hamiltonian 
quantum system in the presen ce of a bath (the exten- 
sion to several baths is trivial) (jBreuer and Petruccionel 

Hoi (n = i), 

vfr(i) = -iWSi{t) - tyvv(t) + l(t)VV(t) . (26) 

Here <J/(i) is the many-body state vector, "H is the Hamil- 
tonian of the system (describing both the molecule, the 




FIG. 11 (Color online) Upper panel: schematic representa- 
tion of the model molecular junction composed of two quasi- 
two-dimensional leads connected with a molecular wire. The 
leads are coupled to external heat baths, each at its own tem- 
perature, (a) Electron charge imbalance as a function of the 
ratio between the couplings gn(L) between the wire and the 
right (left) lead. A strong de pendence can be observed an d 
even a change of sign (from (Du bi and Di Ventral [2009d)). 
(b) Distribution of charge imbalance, when the couplings be- 
tween the wire and the leads are drawn from a Gaussian dis- 
tribution, with an average g — 0.1 and width Ag — 0.001, 0.01 
and 0.1 in units of the hopping parameter (see text). 



leads and the coupling between them) , V are the so-called 
bath operators (which could in principle b e position 
and/or time-dependent) (|Van Kampenl . [200 ll ). which de- 
scribe transitions between the different many-body states 
induced by the bath(s), and l(t) is a stochastic field which 
is taken to have zero mean and a (5-function autocorrela- 
tion, (|(t)) = 0, (l(t)l{t r )) = S(t - t'). 

As a first demonstration, the method was recently 
used by the present authors to study the thermopower of 
a simple model system of spinless non-interacting elec- 
trons, for which the calculations can be equivalently 
carried out with the density-matrix rather than the 
state vector, by averaging over the stochastic realiza- 
tions ( Pershin et all 12008! ). In that model two planar 
leads, each in contact with a thermal bath at a given 
temperature, are connected via a nanoscale wire (see up- 
per panel of Fig. [TTj) . The bath-electron interactions are 
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described by the operators 



L,R 
kk' 



%k?f L ' R (tk)\k){k'\ 



(27) 



where \k) are the single-particle states of the Hamil- 
tonian, f L ' R is the Fermi function containing informa- 
tion on the left (right) bath temperature, and 7^^, 
which are the (inelastic) transition rates between states 
k and k', depend on the bath location (i.e., left or right) 
(jDubi and Di Ventral . l2009df ) . The corresponding equa- 
tions of motion are then solved to obtain the wave- 
function, and hence the electron density, potential, and 
also the heat currents (see also Sec. IIII[) at steady state. 

Several interesting features are revealed by this model. 
For instance, the obtained thermo-voltage shows non- 
linear characteristics, which imply that the linear regime 
may no t be the best regime t o opera te a thermo-electric 
device ( Dubi and Di Ventral . l2009dl) . Another interest- 
ing feature is the strong dependence of the charge imbal- 
ance in the leads (and hence the thermopower) on the 
coupling between the wire (or molecule) and the leads. 
This confirms the ex perimental findings in metallic q uan- 
tum point contacts dLudoph and Ruitenbeekl . 1999). In 



Fig. rjjTa) (|Dubi and Di Ventral . l2009dl ) the charge im- 
balance between the leads is plotted as a function of 
the ratio between the coupling between the wire and the 
left and right leads. From the figure it is obvious that 
the charge imbalance strongly fluctuates, and can even 
change sign as a function of the wire-lead coupling. To 
demonstrate the importance of these fluctuations and to 
tie with the experimental results presented in Sec. IIV.B1 
we have perfor med a calculation fo r the s ame system 
as discussed in ( Dubi and Di Ventral l2009dl) . where the 
coupling constants were drawn from a normal distribu- 
tion around typical values of g = O.lt, where g is the lead- 
wire coupling and t is the tight-binding hopping param- 
eter, which describes the bandwidth of the leads. Other 
numerical parameters are the filling fraction of electrons, 
11 = 1/3, and the temperatures of the left and right heat 
baths, Tl = 0.1 and Tr = 1 (in units of t). A histogram 
of the resulting charge imbalance is plotted in Fig. [TTT b') , 
for three values of the width of the normal distribution, 
Ag = 0.001, 0.01 and 0.1 (in units of the hopping param- 
eter). While more work needs to be done to explain the 
experimental data presented in Fig. IHJc), these theoreti- 
cal results clearly bear some resemblance to experiments 
by showing a double structure in the charge imbalance 
as a function of the coupling asymmetry. 

Using Stochastic TDCDFT one can extend the above 
system model to interacting sy stems, as well as to 
a mu lti-component formulation ( Appel and Di Ventral . 
2009), whereby the Hamiltonian now contains the cor- 
related motion of electrons and (possibly quantum) ions, 
with both components interacting with an external en- 
vironment. Such studies would shed new light on the 
role of interactions and ion dynamics on thermopower, 
and enable a study of local heating effects in nanoscale 
systems (see also Sec. IIII[) . No results are, however, yet 
available for these cases. 



V. SUMMARY AND FUTURE PROSPECTS 

In this review we have discussed energy transport 
in nanoscale systems, such as molecular junctions, sus- 
pended nanotubes, quantum point contacts, etc. Our 
aim was to put under a unified theme the three major is- 
sues of thermal transport, namely thermal conductance, 
local temperature and heating, and thermo-electricity. 
We have critically examined both theoretical and experi- 
mental aspects of these topics. We have presented many 
theoretical methods based on the single-particle scatter- 
ing approach, non-equilibrium Green's functions formal- 
ism, molecular dynamics, etc. From the experimental 
side we have reviewed state-of-the-art experiments, and 
stressed the difficulty and open questions in analyzing 
such experiments. 



A. Future Prospects 

We wish to conclude this review by presenting three 
novel ideas related to energy transport in nanoscale sys- 
tems. These ideas, which deviate somewhat from the 
usual path of thermoelectricity and heat transport, re- 
flect to our opinion the richness and usefulness of study- 
ing energy flow in nanoscale systems, and we hope they 
will stimulate both the experimental and theoretical com- 
munities. 

Thermo- spintronics - Thermo-spintronics (some- 
times also called spin calorimetrics) refers to the ma- 
nipulation of electron spins with thermal effects. Gen- 
erating spin currents, that is the flow of electron spins, 
plays an eminent role in the field of spintro nics, which is 
the s pin analogue of electronics (see, e.g., ( Zutic et~al\ , 
2004)). However, manipulating spins in order to gen- 
erate spin currents is quite difficult, and it is equally 
hard to generate a spin current without generating an 
accompanying charg e curr ent. In recent experiments 
(Uchida et ai, I2008L l2009h a spin-analog to the See- 
beck effect was used to generate a spin-voltage, in- 
duced by a temperature difference along a ferromag- 
netic slab. Although this effect is rather small (com- 
pared to its charge counterpart, but larger than ex- 
pected in view of spin-flip scattering) and inherently in- 
duces an electric voltage as well, it has been suggested 
that these shortcomings may be overcome by applying a 
temperature gradient to a molecular junction placed be - 
tween ferromagnetic leads ( Dubi and Di Ventra . 2009elh 
a setup which wa s further studied recently ( Lii et all , 
\201(i IWang et ali . l20lfl lYing and Jirl |2010D . In an- 
other interesting work, a variety of therm oelectric effects 
in magneti c junctions have been studied ( Hatami et all 
I2007L l2009t iHeikkila et all . l2010l) . with unusual features 
such as thermal spin-transfer torque, spin-polarized cool- 
ing and spin-heat coupling effects. 

Enhanced thermopower in DNA - DNA, the basic 
building block of our genetic code, sh ows also large poten- 
tial in nanotechnology applications ( Dekker and Ratnerl . 
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In a recent study it was shown ([Maria 
20071 ). using a model Hamiltonian for different DNA-likc 
chains, that under certain conditions, the Seebeck coef- 
ficient and figure of merit of a lead-DNA-lead junction 
can be quite high, and can rise to hundreds of fj,V/K, to 
be compared with a few /uV/K of other single-molecule 
junctions studied so far. These high values of the ther- 
mopower seem to stem from transport resonance effects, 
which can be tuned rather easily in DNA. This, with the 
fact that there is a lot of know-how regarding DNA ma- 
nipulation and preparation, makes DNA-based systems 
interesting candidates for future thermoelectric and cool- 
ing devices at the nanoscale. 

Thermoelectricity in superconducting wires - Raising 
the critical temperature, T c , of superconducting mate- 
rials is clearly a technologically important goal. How- 
ever, most superconducting materials have a T c well be- 
low room temperature, even in the well-known high-T c 
materials (with T c ~ 80K for wires and T c ~ 200K 
for bulk). Recently, the pre sent authors have suggested 
(jPubi and Pi Ventral . l2009ah to study a superconducting 
wire held at two different temperatures at its edges. Us- 
ing the method introduced in Sec. IIV.C. 21 combined with 
a self-consistent mean-field theory, it was shown that for 
an (ideally) clean superconducting wire, if one of the sides 
(the cold side) is held at low enough temperatures, the 
temperature of the hot side can be much larger than the 
equilibrium T c , with the wire still in its superconducting 
state. Although this study neglects some effects (such 
as phase fluctuations), the basic idea is simple: in order 
to have a superconducting wire, instead of cooling down 
the entire apparatus, one can locally cool the wire by at- 
taching to it a local refrigera tor, for exampl e one made 
from a Peltier cooling device ( Shakourl 120061) . This may 
pave the road for hybrid superconducting circuits which 
operate at relatively high temperatures. 



B. Final thoughts 

These last few examples - and what we have discussed 
in this Colloquium - clearly show that the quest to un- 
derstand energy transport in nanostructures is far from 
over. In fact, it seems to us that we have barely scratched 
the surface of this problem and more discoveries await 
us. Regarding thermal conductance, finding systems that 
show either very good (for nano-electronic applications) 
or very poor (for thermo-electric applications) thermal 
conductance is needed. As for thermoelectricity, there 
is a need to advance our theoretical tools quite substan- 
tially. For instance, theories that account for the statis- 
tical nature of the experiments should be developed that 
include also electron-electron and electron-phonon inter- 
actions on equal footing. In addition, since the problem 
is intrinsically out of equilibrium (even at steady-state) 
these theories need to include dynamical effects. As for 
local heating and local temperatures, the handful of ex- 



periments that have appeared recently are certainly a 
great start, but more are needed in order to truly de- 
termine the processes leading to heating (and cooling) 
in nanoscale junctions. Similarly, more experiments that 
could determine directly the validity (or invalidity) of 
Fourier's law are highly desirable. 

Pue to the rapid developments in science and technol- 
ogy it is difficult to predict where the field will go from 
here. However, there is no doubt that novel and inge- 
nious ideas will be put forward that will help us profit 
from energy flow, storage and conversion. Embarking in 
such a quest could not be more timely. 
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